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SUMMARY 

A study has been made to determine the feasibility of replacing 
variable-pitch propeller mechanisms with circulation-control propellers . ; 

on general aviation airplanes. The study used a specially-developed 
computer program written in BASIC and placed emphasis on comparing the * 

aerodynamic performance of circulation-control propellers with j 

! 

conventional propellers. ] 

The aerodynamic performance of circulation-control propellers is 
compared with the aerodynamic performance of both variable-pitch and j 

fixed-pitch propellers against the requirements of a 1600 Kg (3600 lbs) j 

| 

single-engine aircraft. The application of a circulation-control j 

propeller with a supercritical airfoil was found feasible under a 

representative design conditions. All propellers had approximately the 
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same performance at high speed cruise (design condition). At low speed, 
the performance of the circulation-control propeller was higher than 
that for a fixed pitch propeller but lower than that for a variable 
pitch propeller. 

It appears feasible to replace variable-pitch propellers with 
circulation-control propellers on single engine aircraft or on 
multi-engine aircraft which have their propellers on a common axis 
(Tractor-Pusher). The economics of these replacements require a study 
for each specific aircraft application. 


1.0 INTRODUCTION 

The feasibility of using circulation-control (C/C) airfoils 
(blunt-based airfoils utilizing the Coanda effect) for general aviation 
airplane propellers Is being appraised. Fuel consumption and/or total 
life cost advantages might occur through replacement of variable-pitch 
propellers with fixed-pitch C/C propellers. The required changes in 
propeller aerodynamic characteristics throughout the speed range of the 
airplane could be obtained through changes In the mass-flow rate of the 
blown jet. A simplified analytical approach with hand calculations 
provided a first-order estimate of propeller aerodynamic performance 
with elliptical and supercritical (S/C) circulation-control (C/C) 
airfoils (ref. 1). These results Indicated that a S/C-C/C airfoil for 
which data became available appeared aerodynamical ly suitable for use In 
a propeller. This study was undertaken to more accurately define the 
characteristics of C/C propellers by a more refined analytical approach 
which used a specially-developed computer program written In BASIC. 


2 


^ & . -> T i ^ 



a* '“ J * 


9b\h 


This program Is presented In the Appendix. A brief assessment was made 
of the required air compressor and the Installation considerations. 

As part of the assessment of a S/C - C/C propeller, performance 
comparisons were made for a selected airplane witf. S/C variable-pitch 
and S/C fixed-pitch propellers. Although none of the propellers were 
optimized for maximum performance, the basic relative characteristics 
should be valid. The propellers were designed for a high-speed 
steady-state flight condition of a typical single-engine airplane. 
Further comparisons were made at a low-speed flight condition for 
steady-state performance and for rate-of-cl Imb or acceleration 
performance. 

This study was initiated by Mr. H. Douglas Garner of the Langley 
Research Center, who proposed the concept as a result of his work with 
fluidic devices. The study received support from Mr. Emanuel Boxer, 
Distinguished Research Associate LaRC and the analytic software was 
adapted from a procedure developed by Mr. William H, Phillips, 
Distinguished Research Associate LaRC, Mr. Wayne H. Bryant of The LaRC 
programmed the technique in BASIC, primarily for this study, however the 
program can support other general design evaluations. The program 
prepared by Mr. Bryant appears as the Appendix. 
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SYMBOLS AND COEFFICIENTS 
Induced axial velocity ratio 
Number of propeller blades 
Propeller chord, m (ft) 
Airfoil wake drag coefficient 
Airfoil lift coefficient 
Blowing Momentum coefficient, 
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m V 


WpIVc 


Propeller diameter, m (ft) 

Propeller aerodynamic horsepower, 746 watts 


550 ft-lbs 
sec 

. . 


Compressor horsepower, 746 watts | 550 ft-lbs 

sec 


550 ft-lbs 
sec 


Propeller-compression horsepower, 746 watts 


T V 

Useful horsepower ” , 746 watts 

55U” 


550 ft-lbs 
sec 


550 ft-lbs' 
sec 


Required horsepower, 746 watts 550 ft-lbs 

sec 


avail Available horsepower, 746 watts 


550 ft-lbs 


sec 


Altitude, km (ft) 

Ratio of Lift Force to Drag Force for an airfoil 
Blowing mass flow per unit span, kg/sec/m (slugs/sec/ft) 
Propeller rotational speed, rev/sec 
Static pressure N/m 2 (lb/ft 2 ) 

Total pressure, N/m 2 (lb/ft 2 ) 

Ambient pressure, N/m 2 (lb/ft 2 ) 


i i 


v 






. ; *ar • - 




ORIGINAL PAGE 13 
OF POOR QUALITY 

Q Propeller torque, N-m (ft-lb) 

R Universal gas constant, m 2 /sec 2 °K (ft 2 /sec 2 °R). 

Also: Propeller radius at the tip of the blade, m (ft) 

(RN) Reynolds number based on chord, C 1 V 1 p » 

P. 

r Propeller radius at a radial station, m (ft) 

T Static temperature, °K (°R) 

Also: Thrust, N (lb) 

V Velocity, m/sec (ft/sec) 

V Free stream velocity, m/sec (ft/ sec) 

W Airplane gross weight, kg (1b) 

a Angle of attack, deg (See Diagram below) 

$ Propeller blade angle, deg (See Diagram below) 

n Propeller Efficiency 

y 

X Advance Ratio, (See Diagram below) 

2irftn , 

P „ Ambient viscosity, N sec/m 2 (lb sec/ft 2 ) • 

p Density, kg/m v (slugs/ft J ) j 

I 

P w Ambient density, kg/m 3 (slugs/ft 3 ) 

♦ Angle of advance, tan “*f ^ °° 1 , deg (See Diagram below) i 

ITV n rl ! 

i 

[ 
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Stibcrl pts 

h Blade hub station 

1 Blade element radial station, (0 to 20) 

(NOTE: Units for T and Q with subscript 1 are per unit span) 

j Jet 


Propeller Blade Definition of Terms 






3.0 RESULTS AND DISCUSSION 


3.1 Choices of Analysis Conditions 

3.1.1 Choice of Airplane 

The design and comparisons of propeller performance were based upon 
applications to a specific model general aviation airplane. A survey of 
general aviation airplanes (ref 2) showed the following ranges for 
maximum speeds. 

(a) Single Engine, Fixed Pitch 59.2-77.1 m/sec (194-253 ft/sec) 

(b) Single Engine, Variable Pitch 66.8-97.6 m/sec (219-320 ft/sec) 

(c) Twin Engine, Variable Pitch 87.5-123.5 m/sec (287-405 ft/sec) 
One of the well established models of single-engine general 

aviation airplanes showed a maximum speed of 82.3 m/sec (270 ft/sec) 
which is about the maximum for which a fixed pitch propeller might be 
used. This airplane was selected because its performance overlapped the 
range for both fixed and variable pitch propellers. The aircraft 
selected had a maximum weight of 1600 kg (3600 lbs) and used a 3 bladed 
variable pitch propeller driven by an unsupercharged engine. Alternate 
models of the airplane are available with supercharged engines. 

Reference 1 presents some of the specific characteristics of this 
airplane. 

3.1.2 Choice of Propeller Airfoil 

A result from the previous study (ref 1) Indicated superior 
performance for a propeller which utilized C/C-supercritlcal airfoils as 
compared with the previously analyzed C/C -elliptical airfoils. 

Available data for S/C - C/C airfoils (ref 3) were therefore used for 
this study. In the previous analyses, preliminary data were provided by 
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the David Taylor Naval Ship Research and Development Center with the 
drag characteristics obtained on a 17 percent chord thick S/C airfoil 
reduced to values equivalent to a 15 percent chord thick S/C alrf' ' 
direct comparison with the previously tested 15 percent chord t: 
elliptical airfoils. For this study, no thickness correction wa* 

The probable use of thinner airfoils in an actual >»opel er application 
Is expected to have an Insignificant effect on the prupeller. 

3.1.3 Choice of Flight Conditions 

A high-speed cruise flight condition of 82.3 m/sec (270 ft/sec) at 
3.05 km (10,000 ft) altitude was selected as the design point for the 
propellers. Although the performance of the airplane should be computed 
at many flight conditions, a selected low-speed flight condition of 38.1 
m/sec (125 ft/sec) at sea level should provide a meaningful Indication 
of off-design propeller performance. Relative propeller efficiency In 
steady-state low-speed cruise and relative thrust margin available for 
climb or acceleration with the given engine provides a useful comparison 
of the propellers. A practical design must evaluate the entire range of 
flight conditions Including take-off performance, obstacle clearance, 
maximum speed etc. This limited comparison Investigated only two flight 
conditions. 

3.1.4 System Description 

The system analyzed. Figure 1, represents the simplest concept that 
could be applied to the circulation-control propeller. Air Inducted 
from an Inlet at ? m Is fed to an engine-driven compressor. This air Is 
furnished to the propeller hub via ducting, valves, and rotary seals. 

The propeller provides additional compression by centrifugal pimping so 
that a variable pressor* head exists from the hub to the outboard 
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Figure 1, Circulation Control Propeller System 
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end of a single blowing-slot plenum chamber. The exit slot Is assumed 
to be the major flow restriction In the system so that the pressure 
gradient In the propeller plenum Is calculated from the hydrostatic 
equations for pressure equilibrium. The radial cross-flows which may 
exist In the exit slot and exhaust have been ignored in this study. 

3.2 Derivations 

3.2.1 Aerodynamic Relationships 

The propeller section aerodynamic characteristics are described by 
the standard c^and cj values as functions of angle of attack a for the 
non-blown propellers. Circulation control propellers require the 
additional parameter of momentum coefficient . These relationships are 
shown graphically In Figure 2 for the supercritical 17 percent thick 
airfoil, and the same airfoil modified for trailing- edge blowing. The 
polars shown were linearized over small angle- of- attack ranges 
(usually 3 degrees) with the numeric results, extrapolated, shown in 
Tables 1, 2, and 3. The software described In the Appendix linearly 
Interpolates between these tabulated values to obtain specific values of 
the aerodynamic characteristic:.. 

The momentum coefficient at any section Is related to the plenum 
total pressure as follows: 


■ v i <«„/ a) 

This relationship Is derived from Figure 3 (reproduced from ref 1). 
Sample points computed from the above equation have been added to the 
original figure. 

10 







Table 1, Characteristics of a Supercritical Circulation 
Control Airfoil (ref 3) 
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Table 2, Interpolation Coefficients for the Characteristics 
of a Supercritical Circulation Control Airfoil (ref 3) 

c t and c d = A + Bc u + (C + Dc u ) f^' ^ 1 


Coefficients to determine c^ 


“1 

lower 

n u 

upper 

A 

8 

C 

D 

-12 

-9 

-0.68 

53.5 

0.26 

-1.0 

- 9 

-6 

-0.42 

52.5 

0.27 

-1.0 

- 6 

-3 

« 

o 

• 

»— » 
cn 

51.5 

0.28 

-2.0 

- 3 

0 

+0.13 

49.5 

0.25 

-0.50 

0 

+3 

0.38 

49.0 

0.28 

-3.50 

+3 

6 

0.66 

45.5 

0.26 

-5.50 

6 

9 

0.92 

40.0 

0.24 

-8.00 

9 

12 

1.16 

32.0 

0.17 

-9.50 

12 

14 

1.33 

22.5 

0.07 

-5.50 

14 

>14 

0 

0 

° 

0 



Coefficients to determine 

c d 


-12 

-9 

0.009 

-0.40 

-0.001 

0.000 

-9 

-6 

0.008 

-0.40 

0.0000 

0.000 

-6 

-3 

0.008 

-0.40 

0.0009 

-0.025 

-3 

0 

0.0089 

-0.425 

0.0011 

+0.020 

0 

+3 

0.0100 

-0. 

0.0031 

0.040 

+3 

6 

0.0131 

-0.365 

0.0024 

0.075 

6 

9 

0.0156 

-0.290 

0.0069 

0.000 

9 

12 

0.0224 

-0.290 

0.0035 

0.020 

12 

14 

0.0259 

-0.270 

0.0077 

0.270 

14 

>14 

sin a 

sin a 

sin a 

sin a 


_ iji 




Table 3 Characteristics of a Supercritical 
17 Percent Thick, Thin-Trailing-Edge 
Airfoil (ref 3) 
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Jet Momentum Coefficient (ref 1) 
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The Btrnoulli incompressible relationship between velocity and 
pressure is: 

P = 1/2P W (V. ) 2 + Po) (2 ) 

This equation gives the jet velocity when the section plenum total 
pressure air is expanded to free-stream static pressure. 

The plenum pressure at any station is related to the hub pressure 
as shown below by the isothermal hydrostatic equation which accounts for 
centrifugal acceleration as: 


dP t • P(2nn) 2 r dr 

(3) 

II 

(4) 


dp _ ( 2 Trn) 2 r dr ( 5 ) 

P RT 


Integrating from the hub to the radius at station i 



( 2 rr n) 2 r 2 


■JlfT 


e 


E i 


( 6 ) 


The hub pressure can be adjusted by use of the control valve which 
operates in series with the compressor. 

The local alrstream velocity at radius r Is the vector sum of the 
forward =md rotational velocities: 

V 1 » [v * 2 + (ZmyO 2 ! (7) 
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This equation does not take into account the inflow velocities 
induced by the propeller but the corrections are too small to be of 
concern. 

Combining the above relationships, an expression can be derived 
between the momentum coefficient and the hub pressure: 


c 


M i 


2RT 

. ] 


1 

(26.4)* 

Poo 

W. m 


v J + (2nr.n) 2 
•» * 


( 8 ) 


The temperature T to be used in this equation is the propeller 
temperature; however, at the low speed conditions evaluated in this 
report, the free-stream static temperature can be substituted. 

3.2.2 Overall Propeller Characteristics 

Propeller thrust and torque are defined by the following 
relationships: 


i = Hr * r P » V » 


1 + a 


sin <j> 


Bc(c£COs<t> - c d sin<t>) 


(9) 


Qi = ir = h'Wj pfinfj Bc(c£Sin<|> + c d cos$) 


( 10 ) 


These are the conventional propeller section equations and account 
for the induced flow at each section. To compute the overall thrust and 
torque, the section contributions are summed using Simpsons rule. The 
method used for computing the inflow by Iteration Is described in the 

Appendix and Is adapted from a method developed by Larrabee and French at 
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the Massachusetts Institute of Technology. 

In addition to the strip integration method described above, use 
was made of an analytic software program developed by William H. 
Phillips, (Distinguished Research Associate, LaRC), to secure first 
app» .imations for efficient propeller chord distributions. These 
relationships are also based ori work by Larrabee (ref 4). In all the 
computations made, the analytic and strip integrations for the same 
propeller agreed to within 1 percent. 

3.2.3 System Efficiency 

The overall efficiency of the propeller system is computed by 
dividing the useful work done by the sum of the work required to rotate 
the propeller and the work required to compress the free-stream air to 
the plenum pressure at each propeller station. 

The required mass flow per unit length of propeller is, (ref 1) 



Vi 


PJV 



(11) 


The horsepower required per foot of propeller radius to pump a mass 
of air from the hub to a radial station is: 


HP 







10.286 



(12) 
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For a propeller with B blades, the expression can be written for 
each radius station r: 


( 



* 10.92 m. B T 



(13) 


To find the overall horsepower required, the section horsepower Is 
summed using Simpsons rule. The horsepower for the engine-driven 
compressor is determined from the requirement that the total mass flow 
Is compressed from the ambient pressure to the pressure at the propeller 
hub. The total horsepower into the propeller consists of the work for 
aerodynamic torque plus the work required to centrifugally pump the air, 
as well as that performed by the compressor to supply hub pressure. 


HP 


total 


= HP + 
aero 


HP 


pc 


HP. 


HP total 


2imQ 

550 


+ HP 


pc 


+ HP„ 
c 


(14) 


The useful work done by the propeller is the thrust multiplied by 
airplane velocity, and In horsepower Is: 



TV. 


(15) 
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The system efficiency is: 
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TVqq/550 

SP* H V * HP c 
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The above equations, and other ancillary relationships to define 
atmospheric properties versus altitude, local Reynolds and Hach numbers 
and other pertinent factors have been mechanized in BASIC and are 
described in the Appendix. The program is useful for designing any 
propeller for which the section aerodynamic characteristics can be 
stated versus apgle of attack in a look-up table. In addition, the 
program is capable of varying parameters such as engine rotational speed 
or blade helical angle and provides cross-plots for trend analysis. 
Pertinent results of the computations made for this study are presented 
below. 

3.3 Performance 

The design point for all of the propellers was 82.3m/sec (270 
ft/ sec) at 3.05 km (10,000 ft) altitude and they were then analyzed for 
low speed flight at sea level. No attempt was made to optimize the 
propeller diameter or number of blades. A 1.83 m (6 ft) diameter 3 
bladed propeller was used for all comparisons. The listing below 
summarizes the design and off-design flight conditions chosen plus 
pertinent airplane and engine parameters: 
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Alt1tude,h 

Velocity, V w 

Thrust required, 
steady state 

Engine rpm, 

(full throttle) 

Engine HP, 

(full throttle) 


Design (Cruise) 
3.05km (10,000 ft) 
82.3 m/sec (270 ft/sec) 
1441N (324 lbs) 

2500 

188 


Off -Design (Sea Level) 
0 

38.1 m/sec (125 ft/sec) 
1544N (347 lbs) 

2700 

285 


3.3.1 Airfoil Characteristics 

The aerodynamic characteristics measured for S/C-C/C airfoil 
(configuration 5 of ref 3) are shown In Figure 2 for a range of blowing 
momentum coefficient and angle-of-attack. The drag coefficients 
presented are wake drag coefficients that Include the momentum of the 
blown jet. These differ from the coefficients used In the Initial study 
(ref 1), where, for purposes of airfoil comparison, the wake drag 
coefficients were converted to drag coefficients that Included a drag 
equivalent of the blowing power required. In this study, the power 
required to blow the air Is separately accounted for and charged to the 
airplane engine. This method takes account of the trigonometric 
relationships required to compute propeller efficiency. 

In the angle-of-attack range of high L/D (about 3° for the S/C-C/C 
airfoil), the L/D at no blowing Is 53.8 and Increases with blowing. 

The overall efficiency of the airfoil Is limited, however, by the energy 
required to compress the air from ambient pressure to the plenum 
pressure at the jet. As a typical example of section efficiency, the 
characteristics of the S/C-C/C section were evaluated at an advance 
ratio of 0.69 for a fixed 3° angle of attack with various values of Cy . 
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The section efficiency at Cy = 0 of 0.96 was reduced to 0.95 as the Cy 
Increased from 0 to 0.02. Larger values of CyWill further decrease 
the section efficiency even though the aerodynamic L/D for the section 
itself is increasing. For applications where the jet velocity is not 
otherwise limited, it is possible to use momentum coefficients up to 
about 0.02 without excessive losses in efficiency. 

3.3.2 Design Methods and Design Factors Considered 

The number of blades and propeller diameter were held constant to 
study the relative performance of the propellers described below. The 
diameter was selected to maintain a low tip speed which permitted an 
additional jet velocity near the tip. In each case, an angle-of -attack 
distribution was selected and the analytic minimum-loss chord 
distribution determined. The chord distribution was scaled to achieve 
the required thrust at the high-speed design point. 

Three angle-of-attack distributions were investigated: a nominal 2° 

and 4° , uniform hub to tip, and a "twisted" distribution of -12° hub to 
+5° at the tip. The small angles of attack were selected because they 
were located in the region of maximum section L/D and provided for a 
large range of operation before stall at low forward speed. The 
"twisted" distribution was selected to give reasonable performance at 
the high-speed design point and to minimize stall as much as possible at 
the off-design point. The summary below outlines the design procedure' 
used to match propeller performance to the airplane thrust requirements. 








, W 
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Propeller Type 


Design Point 
3.05 km (10000 ft) 


Off-Design 
Sea Level Cruise 


S/C 


Analytic minimum-loss Fixed Pitch - vary 

chord distribution engine rpm 

Variable Pltch- 
vary hub setting 


S/C, C/C, with 
blowing plenum 
to 0.7R and no jet 
velocity restriction 


Match c* at 0.7R to Set Cy * 0. Vary 
above S/C design. engine rpm 
Ratio chords to 
secure required thrust. 


S/C, C/C with blowing Analytic minimum Set Cy ■ 0. Vary 

plenum to tip and loss chord engine rpm 

0.9 sonic velocity distribution 
restriction 


Each case In the parametric study (10 total) produced a print out 
of the propeller characteristics, plus jet characteristics for blown 
propellers. Table 4 shows a typical printout. In addition, the program 
generated a number of cross plots. Figure 4 shows examples. The 
Appendix contains print-outs of propeller characteristics for the 
various designs studied. 

For the off-design condition, the propeller-engine combination was 
evaluated to determine the excess thrust which could be made available 
above that required for steady-state flight. The excess thrust, useful 
for rate of climb or acceleration, was limited by the full-throttle 
engine horsepower at 2700 rpm or by the maximum thrust capability of the 
specific propeller at any rpm up to 2700. 

3.3.3 Propeller Characteristics 

A summary of the operating characteristics of the final C/C and 
ui. i propellers (VP and FP) Investigated Is presented in Figure 5 and 
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* Evaluation of a Fixed Pitch Propeller ** Evaluation of a Variable 

at various engine speeds at various blade settings 

Figure 4. Computer Generated data for Comparison Evaluation of Propellers 
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the propeller efficiencies In steady-state cruise are presented In 
Table 5 for both the high-speed and low-speed flight conditions 
selected. 


Table 5. Calculations of Propeller Efficiency 


Propeller Type 


Variable Pitch 
Narrow Chord 


Design Point 
82.3 m/sec at 3.05 km 
(270 ft/sec at 10,000ft) 


Off Design Point 
38.1 m/sec at Sea Level 
(125 ft/sec at Sea Level) 


0.897 at 2500rpm 0.784 at 2700rpm 


Fixed Pitch 0.897 at 2500rpm 0.620 at 2100rpm 

Narrow Chord 


S/C-C/C 0.886 at 2500rpm 0.777 at 1800rpm 

Wide Chord with 
blowing to tip, subsonic 
jet 





* 


For the high-speed cruise design condition (Figure 5a), the chord of 
each propeller (C/C and unblown) was established as described above 
(section 3.3.2) to provide the thrust required for cruise at the 
recommended engine rpm of 2500 for sustained steady-state flight, 
(Indicated as®). A value of c jj s 0.0025 at 0.7 R was selected for the 
C/C propeller prior to establishment of the C/C chord to maintain 
subsonic blowing. The printout as Table 4 shows the variation of c p 
(and jet velocity when expanded to ambient pressure) along the span. 

The efficiencies at high-speed cruise are nearly the same (about 0.89) 
for all three propellers (see Table 5). This Is not unexpected as all 
are operating at a selected angle of attack that provides high L/D 
ratios. 
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Thrust , N (lbs.) Thrust , N (lbs. 


1778 

(400) 


889 

( 200 ) 


All Propellers 

— 


Thrust Required''^ y ^ 1 
W-l 600kg (3600 lbs) x P 

z' x ./l 
s' Thrust Required 
W«1 333kg (3000 lbs) 


ORIGINAL PAGE IS 
OF POOR QUALITY 

n = 0.39 


S/C-C/C Wide Chord 
© - 0.0025 at 0.7R 

0c„. o 


(a) High Speed Cruise; 82.3 m/sec at 3.05 km, (270 ft/sec at 10,000 ft) 
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(b) Low Speed Cruise: 38.1 m/sec. at Sea level, (125 ft/sec) 


Figure 5, Propeller Operating Characteristics 
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With this propeller efficiency assumed, the thrust available with 
the engine operating at full throttle was calculated and plotted against 
rpm as the dashed line In Figure 5a. Each propeller is capable of 
providing the thrust required at about the maximum horsepower available 
from the selected airplane engine at cruise rpm. Blowing control can be 
used for steady-state high-speed flight at a constant engine rpm to 
accommodate some decrease In airplane weight below the design weight 
(shown asQln Figure 5a, to W =1333kg, 3000 lbs). A further reduction 
in weight, of course, will require reduced engine rpm. In contrast, a 
variable-pitch propeller can opeiate through a wide range of weight at a 
constant engine speed. 

For the selected low- speed condition, the thrust available with 
the S/C - C/C propeller designed with a subsonic jet at the high-speed 
condition Is plotted against rpm in Figure 5b as the solid line. The 
propeller efficiency at the rpm at which the required steady-state 
thrust Is obtained (1800 rpm) and the propeller efficiency at the 
maximum rpm of 2700 rpm are Indicated at the ends of the curve as 0.78 
and 0.55, respectively. The dashed lines show the thrust available with 
the engine operating at full throttle plotted against rpm for various 
assumed values of propeller efficiency from 0.65 to 0.45. 

It Is first seen that more than sufficient engine power Is 
available for steady-state cruise at the thrust-required rpm of 1800; 
l.e., the S/C - C/C propeller efficiency Is greater than that necessary 
(<0.55) for the engine at full throttle to provide the required thrust 
for steady-state flight. The engine, therefore, would be operated at 
reduced throttle for steady-state cruise at this low speed. A margin Is 
then available for climb or acceleration. 
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The maximum rate of climb at constant speed Is directly 
proportional to the difference In thrust available at the maximum engine 
power or the maximum thrust available from the propeller from that 
required for steady-state flight. From Figure 5b, it Is seen that the 
S/C - C/C propeller efficiency is high enough (0.55) at the maximum rpm 
of 2700 to absorb the maximum engine power at that rpm. The thrust 
margin for climb or acceleration which results is 1779N (400 lbs). 

The thrust produced with the S/C propeller operating as a variable 
pitch propeller at 2700 rpm for various changes in hub blade angle, 
are indicated "x" ana the corresponding propeller efficiencies are 
listed. The maximum thrust available with a fixed-pitch propeller 
( AB = o) is less than that obtainable with the S/C -C/C propeller. The 
thrust is limited by blade stall and not by engine power availability as 
indicated by comparison of the FP propeller efficiency (0.52) with the 
full -throttle thrust-available curves. The VP propeller with a hub 
blade angle change of -5° produces the largest excess thrust. For this 
case also, the engine would be operated at part throttle, whereas the 
S/C - C/C propeller operates at full throttle because of a lower 
efficiency (0.55 vs 0.66). 

The results shown in Table 5 indicate that, to match the required 
thrust for steady-state cruise with a fixed pitch narrow chord propeller 
at the selected low speed, an engine rpm of 2100 is required. Although 
the analysis was done at an rpm of 2700 for the VP propeller, it is 
clear from Figure 5b that the VP propeller can be adjusted to operate at 
any preferred rpm. 

Also plotted on Figure 5b are results for a S/C - C/C propeller 
with a narrow chord similar to the S/C FP and VP propeller chords. The 
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C/C propeller with narrow chord wns unable to achieve the required 
low-speed steady-state flight thrust at any engine speed because of 
blade stall at angles of attack less than that for stall of the thin 
trailing edge S/C airfoil. Even If some change In design could permit 
achievement of the required thrust (at about 2100 rpm), no excess thrust 
for climb or acceleration would exist. 

3.4 Operating Boundaries 

The aerodynamic data used (reference 3) were measured over a 
limited Reynolds number and Mach number range. For this study, the c^.c^ 
and Cy vs a relationships were assumed to be invariant. Extrapolation 
of the experimental data over a range of negative angle of attack was 
required for a portion of the study. Operating boundaries that should 
be observed are illustrated In Figure 6 and discussed as follows. 

3.4.1 Use of Positive Lift Coefficient 

It is possible to operate over a small range of negative c^for the 
S/C - C/C airfoils and secure positive c^ by blowing. This should be 
avoided to prevent loss of propeller effectiveness In event the blowing 
Is disrupted by equipment malfunction. 

3.4.2 Use of Low Angle of Attack 

On the blown airfoil, the degree of control of c^with changes in 
the momentum coefficient varies widely with angle of attack. For 
example, at zero degrees angle-of -attack, a change of Cyfrom 0 to .01 
produces a change In cyfrom 0.38 to 0.90, a ratio of 2.37:1 Increase. 

At an angle of attack of +9 , for the same c y change, the c^ varies from 
1.16 to 1.48, a ratio of 1.28:1. This Is almost a 2:1 change In 
effective control by blowing. Because of the desire to work at low 
angles, It was necessary to provide wide chords to provide the required 
28 
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Figure 6, Drag Polars for S/C and S/C-c/c Airfoils (configuration 5 of ref 3) 
with Propeller Operating Boundaries 
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thrusts. As shown above (3.3.3,) only the wide chord C/C propellers 
provide reasonable performance. 

3.4.3 Avoidance of Supersonic Jet Velocity 

For the airplane and propellers studied, propeller tip velocities 
were subsonic, and it was possible to design propellers with circulation 
control out to the tip. Although some computations were made which 
resulted in supersonic jet velocities (e.g., results presented in Figure 
5b), there are no data regarding the aerodynamic performance with high 
jet velocities. To keep the jet velocities subsonic, the momentum 
coefficients had to range from no more than 0.002 at the tip to 0.004 
Inboard. The approximate boundary for sonic momemtum coefficient is 
shown in Figure 6 far the single-plenum propeller. 

3.4.4 Section Efficiency 

The c^/cj plot indicates that maximum L/D is secured for the S/C 
airfoil at angle of attack of about 5°, and for the S/C - C/C at about 
4°. The L/D does not linearly affect propeller performance as It does 
for wings on aircraft. In fact, until the L/D decreases to about 30, 
there are only small effects on propeller performance. The L/D for the 
S/C-C/C Is above 30 for angles-of-attack down to zero degrees. Lower 
angles will have significant adverse effects on propeller performance. 

Momentum coefficients below 0.02 are desired, (sec. 3. 3.1), to 
prevent excessive deterioration of section efficiency. In general, the 
supersonic jet boundary restricts operation to below c y * 0.005} 
however, with multi-plenum systems, it should be possible to operate 
with somewhat higher momentum coefficients than 0.005 near the hub as a 
means to Improve propeller performance. 
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3.5 System Design Considerations 

In this limited study, only aerodynamic performance was computed 
and analyzed. Other design considerations such as noise, structural 
analysis and design, system weight, ducting and seal design, and type of 
compressor were considered but not in enough depth to warrant many 
conclusions; however, no large technical problems were apparent. 

Because of the small slot size, it is judged that the noise energy 
would be located at high frequency and therefore not be a problem. 

Also, blowing is only used for the high-altitude high-speed condition 
and not at low altitudes. The structural problems were only examined to 
the extent that it appears reasonable to make the trailing edge assembly 
as a sub-assembly of the main propeller blade. The small size and 
dimensions of the slot lead to concerns regarding nicks, distortion and 
structural failure under concentrated loads. The design of pneumatic 
shaft seals is believed to be straightforward. 

The compressor requirements to compress the ambient air to the hub 
pressures required for a typical propeller are listed on Table 6. While 
the compression ratio is higher than that required for supercharging the 
airplane engine, the mass flow is small and the compression horsepower 
only about one fourth of that required for the engine alone. Valving of 
the charger outputs would be suitable to control the hub pressure. If 
the airplane engine is not supercharged, it is possible to operate an 
automotive type positive displacement unit (vane or Roots type) as a 
belt-driven accessory. The cruise intake air pumping requirement is 
0.92 liters (56 in 3 ) per crankshaft revolution and is comparable to the 
pumping volume rates available from the automotive market. 
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Table 6. Air Flow Requirements For a S/C-C/C 
Propeller Driven by a Supercharged Engine 



Propeller* 

Engine 

Total 

Air Flow Required 

0.038 kg/sec 0.0159 kg/sec . 0.097 kg/sec 
(0.0026 sl/sec) (0.0310 sl/sec) (0.0336 sl/sec) 

Compression Ratio 

1.22 

1.078 

- 

Compressor Horse- 
power Required 

1.214 

5.02 

6.23 

Compressor Horse- 
power Available 



18.2 


Note: Operating conditions for 3.06 km, (10,000 ft), at 2500 rpm with 
a Turbocharged engine of 8.52 liters, (520 in-*), 
displacement, (ref 2) 


♦The propeller requirements include an allowance for leakage and 
losses within the ducting. 



4.0 CONCLUDING REMARKS 


A specially-developed computer program (presented in the Appendix) 
has been used to compare the aerodynamic performance of 
circulation-control (C/C) propellers with variable-pitch and fixed-pitch 
propellers. The comparisons were made for a 1600 kg (3600 lb) 
single-engine general aviation airplane with a maximum speed of about 
300 km/hr (186 miles/hr), the approximate limit for a fixed-pitch 
propeller. 

The study indicated that, on an aerodynamic performance basis, the 
circulation-control propeller is feasible. Increased speed decreases > 
the potential feasibility. Inability to feather and reverse thrust 
limits applicability to single- engine airplanes or multi-engine 
configurations where engine failures can not produce a disturbing torque 
(e.g. on-axis configurations). Economic feasibility requires analysis 
of manufacturing and maintenance costs of C/C propellers as well as 
appraisal of mission requirements for specific airplane applications. 

All of the propellers investigated had approximately the same 
efficiency at the high-speed cruise design condition. At low-speed, the 
C/C propeller performance (cruise, rate of climb, and acceleration) was 
better than that of an unblown fixed-pitch propeller but not as good as 
that of a variable-pitch propeller. Although blowing at high-speed 
permits operation through a wider range of angle of attack than for an 
unblown fixed-pitch propeller, performance is constrained by the amount 
of blowing permissible. The amount of blowing is limited by decreases 
in efficiency with increases in blowing power and by the desire to 
maintain the blown jet velocity to subsonic values. The latter 


33 









4 



i 


constraint was Imposed in the Interest of conservatism because no 
experimental data are available with supersonic jet velocities. It 
appeared reasonable to expect that the aerodynamic effectiveness of the 
Coanda jet would deteriorate with supersonic blowing. Improved 
performance, however, appears possible through compartmentation of the 
blowing plenum along the propeller span to provide increased subsonic 
blowing at the lower-speed inboard sections. The flexibility of 
operation for a variable-pitch propeller will likely yield superior 
performance at all off-design low-speed conditions. The possible 
overall advantages of C/C propellers, therefore, depend upon economic 
comparisons. 

A limited appraisal of other than aerodynamic design 
considerations, such as noise, structure, weight, ducting, and seals, 
indicated no large technical problems for C/C propellers. The 
compressor requirements can be met with automotive-type compressors or 
by a small amount of valve-controlled bleed from an engine supercharger. 
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APPENDIX 

by 

Wayne H. Bryant 
NASA Langley Research Center 

Introduction 

This appendix describes the computer program developed to generate the 
data in the present study. This aerodynamic propeller design program can 
accommodate minimum loss (ref. 4) fixed- or variable-pitch propellers, and 
circulation-controlled propellers with a single plenum extending from the root 
to any given radius. The program was implemented in BASIC so that only a 
minimal computer investment is required to use it to design propellers. The 
present work was accomplished using a Digital Equipment Corporation VAX-11/780 
superminicomputer, and the BASIC language elements conform to those used in VAX 
BASIC. An effort was made to develop as much of the program as possible in 
"standard" BASIC to ease its transition to other machines. To test this, the 
program was transferred to a CDC Cyber 175 computer at LaRC where approximately 
two hours were required to obtain successful operation. Comments on probable 
coding changes required/de sired appear at the end of this appendix. 

The results obtained using the computer program described here have been 
examined and appear reasonable. While there are no known problems or "bugs" 
remaining in this program, there may yet be problems that will surface for 
new input cases. Additionally, the program makes no structural analysis of 
the designed propellers; the structural integrity must be ascertained by the 
propeller builder using some other technique. 

This appendix is organized into seven main sections. These ares 
Introduction 

General Description of Major Program Sections 

Case Table Identification for the Present Study 

Program Variable Definitions and Program Listing 

Detailed Program Description (keyed to listing line numbers) 

Suggestions for Tailoring Program 

Concluding Remarks 






36 



General Description of Major Program Sections 

The purpose here is to give a brief description of each major area 
encountered in the program. This will be done first in the order the sections 
are found in the listing itself without regard to the program execution flow. 
Next, a specific case will be given to illustrate typical program flow, and 
appears in the Detailed Program Description section. 

There are eight main program sections: 

(1) Case selection 

(2) Minimum induced loss propeller design (analytic) 
and output 

(3) Strip integration propeller design (Part 1), 
both non-blown and blown (circulation controlled) 

(4) Induced velocity iterative calculations 

(5) Strip integration propeller design (completion) 
including thrust matching iterative calculations 

(6) Miscellaneous calculations (e.g ., efficiency, Mach number, etc.) 

(7) Line printer and cross plot file output 

(8) Subroutines 

The main purpose of each section is given in the text that follows. The line 
numbers listed correspond to the program listing found later in this appendix. 
While this program was developed to generate data for the study presented in 
the main body of this paper, it should be remembered that relatively simple 
modifications to the listed program render it useful for evaluating a variety 
of propeller designs. The data contained in the program, and the logic state- 
ments controlling the program flow, are those used for the last part of the 
study. 

(1) Case Selection 

After the array declarations and opening output files, the first part of 
the program consists of initializing variables for a specific run. This occurs 
between lines 170 and 1820. A number of DATA statements contain information 
for the blown propeller lift and drag coefficient lookup table as well as that 
for ten pre-defined evaluation cases. These DATA entries are identified in the 
Detailed Program Description section following the program listing and cross 
Reference table. Table A1 in the next part of this appendix shows the main 
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characteristics of these pre-defined cases. If a run is desired for which no 
case exits, all required data can be manually entered from the keyboard during 
program execution. All propeller designs usi *, manually entered data will be 
evaluated as "on design point" cases; that iu, the propeller's performance in 
some arbitrary off-deeign point cannot be made. The "off-design point" eval- 
uation feature is built in with the use of pre-defined cases. 

The engineering units for each input parameter is displayed ar the value 
is requested. These seme units are used for pre-loading the defined cases in 
the data statements. The data order for the defined cases will be given in 
more detail later. In addition to data for the actual propeller design, the 
program inquires whether either of two forms of diagonistic output is required 
for the current run (both predefined or manual input cases). The first set 
of output is routed to the printer and is useful for observing the convergence 
during induced velocity and thrust matching iterations; the second set routes 
the inflow iteration data to the console display device (CRT) so that the user 
can observe convergence in real-time. 

(2) Minimum induced loss propeller design and output (Analytic) 

The next section designs a non-blown propeller at the given operating 
condition using the technique described by Larrabee in reference 4. This 
section is based almost entirely on a propeller design program developed by 
H. Hevitt Phillips and E. Eugene Larrabee for an HP 9830 desktop computer. 

The code implemented ranges from line 1820 to line 3860; equations noted 
in the REMark statements refer to the numbered equations of (ref. 4), and are 
all contained in the listing within angle brackets (e.g. < EQN 21>). 

After loading data required for the non-blown propeller table look-up 
subroutine, this section obtains the lift and drag coefficients for the angle 
of attack spanwise distribution previously entered. This is accomplished 
through a subroutine call to the non-blown propeller lift/drag subroutine. 
After these quantities are calculated, this section determines the spanwise 
chord distribution, efficiency, required horsepower and torque, the local 
Mach and Reynolds numbers, and other pert inant data. These data are then 
output to the line printer file. 

At the end of this section, chord and beta spanwise distributions are 
saved for later predefined case analysis. A more complete description of this 
mechanism is given in the detailed program description. 
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(3) Strip integration propeller design (Part 1) 

This section is coded in lines 3870 through 4590 and implements equations 
fcuud in the main body of this report. These references are denoted within the 
REMark statements using square brackets (e.g. [EQN. 7] at line 3900). The 
strip integration section can design either blown or non-blown minimum induced- 
loss propellers; this feature is controlled by the designer at program setup 
time. 

The strip integration designs were checked against the analytic designs 
and gave results identical within a few percent although there were small 
changes in the spanwise chord and angle of attack distribution. These changes 
in angle of attack result in an increased thrust which is compensated for (in 
a later program section) by scaling the spanwise chords to achieve the 
required thrust. Comparison of the output from the (non-blown) analytic design 
with that of the (non-blown) strip integration design shows small changes in 
the spanwise chord and angle of attack distributions. 

(4) Induced velocity iteration calculations 

This section implements the iterative equations necessary to determine the 
true spanwise blade angle of attack which is different from the nominal because 
of induced velocity. Lines 4600 through 5420 realize this procedure. The 
technique used here is derived from (Ref. 5) and is listed here in full for 
convenience. References to equations in the comment portion of the individual 
statements or in REMark statements are enclosed in square brackets and have an 
"A” (for appendix) prefix. For example, line 4710 has a comment (the text 
following the "!") indicating that line implements equation A1 ([EQN. Al]); 
that equation itself follows. 

The induced velocity components are evaluated at each radial blade station, i, 
by the iterative procedure listed below. To start the iterations 

♦i - tanl (~r-) 

is used to find the initial estimate of as 

a i ' h ' *t (A2) 

The procedure continues by finding the multiplier 
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ORIGINAL PAGE FJ 
OF POOR QUALITY. 



The induced velocity components and are next used to calculate an 
updated $ as 



which is used in obtaining a corrected a in 

a = a- + ) CA7) 

c i 1 1 c i 

The corrected angle-of -attack is used to calculate a corrected $> as 

= Pj * a c ( AR ) 

i i 

which is used if further iterations arc necessary. 

The convergence tests consists of comparing the current iteration value of 
<* c with the previous iteration value of at each blade station and declaring a 
converged solution if the absolute difference at ever)- blade station is less than 
that of 0.05°. Stated mathematically, the test is 

(Oi " « c ) < 0.05° For all blade stations, i. (A9) 
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If any station fails to meet this criteria, the procedure iterates start- 
ing with Equation A3; if the corrected angle-of-attack at all blade stations 
pass the test, program flow is to the second part of the strip integration 
design. In either case, new lift and drag coefficients, based on the most 
recently corrected angle-of-attack, are calculated for both blown and non-blown 
designs as appropriate. 

(5) Strip Integration propeller design (completion), and thrust matching 

This section is coded between lines 5430 and 6520, with the thrust match- 
ing part between 5660 and 6510. After the initial induced velocity corrections 
are made to each blade station's angle of attack, the differential thrust and 
torque are determined. The total thrust is obtained by numerical integration 
and is compared to the required thrust for the indicated flight condition. If 
the developed thrust does not match the required thrust within 1Z, then some 
form of adjustment is employed. This is accomplished in an iterative fashion, 
with the adjustment parameter selected by a combination of case identification 
and control variables. After the adjustment has been made, program flow 
returns to the induced velocity iterations to account for the present thrust 
matching changes. This sequence is repeated until the developed thrust is 
within limits indicated above, or the designer manually terminates the run. 

This program section largly controls the manner in which cases are inter- 
preted, i.e., the sequence in which cases must be run to obtain meaningful 
results. (This sequence is discussed in the detailed examination of the list- 
ing.) For the current program, alpha, beta, engine speed, or the propeller 
chord is scaled obtaining desired thrust according to case under evaluation. 

To modify the program for one's particular needs, a good understanding of this 
section is required, since most changes will be made here. 

(6) Miscellaneous calculations 

This section starts at line 6530 and ends at line 7290. Differential 
torque is numerically integrated to obtain total torque. If the designed 
propeller uses Coanda blowing, then the total air mass flow, and the 
centrifugal and compressor engine horsepower for this air flow are obtained. 
Next, aerodynamic torque is converted to horsepower using 

Qaero " TT§ (A10) 
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! ‘HPI 


2 


Pa jx.,. 


and the total useful work (in horsepower) is calculated. From this data, 
propeller efficiency is then calculated. Finally, at each blade station; 
(*) local Mach number as 


w. - 

with V g 


'i 

(All) 

* velocity of sound at current altitude. 


(b) Reynold's number 
(RN). - p °° C i V i 


with p * air density and y = viscosity 


(A12) 


and (c) the Drag/Lift ratio 



(A13) 


are calculated. 

The section ends with the determination of the available horsepower at the 
current engine speed for the aircraft selected in this study, and the possible 
rate-of-climb for the propeller/engine/aircraft system. 

(7) Line printer and cross-plot file output 

The next program section runs from line 7300 to line 8400. If the case 
under examination calls for off-design point evaluation, the program saves 
pertinent information in a disk file for later processing. The saved data 
and the file format are discussed in the detailed program description. Next, 
data for the appropriate design technique is output. Three forms of printer 
output are available: (a) analytic (non-blown), (b) strip integration (non- 

blown) , and (c) strip integration (blown). The particular combination of 
output printed is determined by case number and control variables. 

After the output is complete, the final beta and chord values for the 
design case propellers are saved for use in the off-design case evaluation. 

This section ends with a query to either examine another case or stop the 
program. 
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(8) Subroutines 

/ 

The program subroutines can be found starting at line number 8410 and 
running to the end of the program at line 9710. Aerodynamic coefficients 
/. for the non-blown propeller are obtained from a lookup table subroutine 

between lines 8410 and 8770. A Simpson's Rule integration scheme is found 
£ starting at line 8780 and ending at line 8920. An atmospheric characteristics 

■'[ _ subroutine, reproducing the values given in the 1962 NASA standard atmosphere 

* report (Ref. 6), is next, running from lines 8930 to 9190. 

The Coanda effect, blown trailing edge propeller aerodynamic coefficients 
subroutine is located starting at line 9200 and ending at line 9460. The 
subroutine between 9470 and 9630 is used to calculate rank intervals for this 
lookup table. The final subroutine runs from 9640 to 9710 and calculates 
the available engine horsepower as a function of engine speed at full throttle 
for the engine used in this study. 

Case Table Identification for the Present Study 

The current program is structured to design six propellers in cases 1, 2, 
■ 3, and 4. These propellers are denoted A, A', B, C, C', and D. Cases 5-10 

are for off-design point evaluation. Table A1 lists the characteristics and 
control parameters for each case. Definitions of the variables and control 
parameters listed at the bottom of each column, and in the text that follows 
; the table can be found in the list of symbols that preceeds the program 

listing. 

Table Al. Propeller Design Program Case Identification 
Evaluate (?) 

Case # Alpha Design Blown/Non-blown 


1 (A/A') +4 

On 

Y(A)/Y(A' ) 

2 (B) 

-12,5 

On 

Y(B)/N 

3 (C/C') +2 

On 

Y(C)/Y(C* ) 

4 (D) 

-12,5 

On 

Y(D)/N 

5 (A) 

+4 

Off 

Y(A)/N 

6 (B) 

-12,5 

Off 

Y(B)/N 

7 (C) 

+2 

Off 

Y(C)/N 

8 (D) 

-12,5 

Off 

Y(D)/N 

9 (A') 

+4 

Off 

N/Y(A' ) 

10 (C) 

+2 

Off 

N/Y(C' ) 

T8 

A(1,I) 

* 

C7 


Fixed During Thrust Correction (?) 

Blowing Chord Beta 
On Y/N Y 

On Y /- Y 

On Y/N Y 

On Y/- Y 

Off Y /- Y 

Off Y/- Y 

Off Y/- Y 

Off Y/- Y 

-/Y N 

-/Y N 

* C5 B7 











Parameter variation by case number identification is controlled by the 
variable listed at the bottom of each column. For those columns with *, the 
parameter is embeded within the program itself, and cannot be changed without 
modifying the source code. The propeller identification for the ten cases set 
up in the accompanying listing are shown under the "Case if” column in 
parenthesis (e.£>. (A/A')). 

An example is helpful in understanding this table. For Case #1, (T8=l) 
two propellers are designed, denoted A and A' . Both propellers have a nominal 
value of alpha of +4 degrees (A(1,I)), and are considered to be designed for 
the operating conditions specified in the Q(l,*) array as indicated by the "On" 
in the Design column. Both a blown and a non-blown propeller will be evaluated 
(C7=3), and for the blown propeller, the trailing edge jets will be active 
(Blowing?). The propellers are designed to match the available thrust to the 
required thrust specified in the Q(l,3) element. This can be done several 
ways. The next four columns specify for each propeller/case which scheme is 
used to match the available thrust to the required thrust. For the blown 
propeller, the alpha valu.s are scaled to accomplish the matching (as can be 
seen by the N(o) in the first part of the entry under "Alpha".) For the non- 
blown propeller, the chords are scaled, as can be seen from the /N (o) under the 
"Chord" column. The variables C5, B7, and C6 are used to set up the Chord, 
Beta, and Alpha columns, and program logic is used to setup the RPM column. 

The predefined cases must be run in a specific sequence to make certain 
that necessary data is available for each case. For the current program logic, 
the sequences are as follow: 

Case #1 Case #3 


/ \ / \ 


Case if2 Case if 9 

Case #4 

| 

Case if 10 

j 

Case if 5 

| 

1 

Case if! 


Case #6 

| 

Case if 8 


Case #1 designs two propellers, A' 

(non-blown) 

and A (blown) at the 


design point. A* is designed first, and the chords are scaled to match 
the required thrust. A is then designed using the final chords obtained 
from propeller A’j thrust matching is achieved by changing the nominal angle 
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of attack. Case #2 designs one propeller, B (blown). The chords used in 
this design are those obtained from the A' propeller, and thrust matching 
is carried out by changing the nominal angle of attack. 

In a similar manner, Case #3 designs two propellers, C* (non-blown) and 
C (blown), also at the design point. Again, C' is designed first and the 
final chords obtained during the thrust matching is used at the design chord 
values for the blown propeller (C) with nominal angle of attack modified to 
obtain the required thrust. Case #4 designs the sixth propeller, D (blown) 
which uses the chords from the C’ propeller, and also changes the nominal 
angle of attack to secure the required thrust. 

Cases #5, 6, 7, and 8 evaluate the four blown propeller designs (A, B, 
C, and D) at an off-design point. Tor these cases, the blowing is shut off, 
and thrust matching is achieved by changing engine speed. Cases #9 and 10 
evaluate the two non-blown propeller designs (A' and C') at some off-design 
condition, and vary the propeller pitch to secure the required thrust. 

Engineering and Program Variable Definition, and Program Listing 


ENGINEERING 


total 


PROGRAM 

A(1,I) 

A(2,I) 

A(3,I) 

A(4,I) 

A$ 


USAGE 

Alpha at each blade station, degrees 
Non-blown propeller lookup routine interval 
constants 

Non-blown propeller lookup routine alpha 
increment, degrees 
Not used 

Output format string 

Alpha at each blade station during manual 
setup of spanwise distribution, degrees 
beta increment/iteration during off-design 
analysis, degrees 

Total beta change during off-design analysis, 
degrees 

Logical: -1 = forced case end during off- 

design analysis 

Square of jft velocity (temporary variable) 
Alpha increment/iteration during off-design 
analysis, degrees 


•'W .i A’iSC-tV 









k... \ 


ENGINEERING 

Act , 
total 


B 


% 

T (Ref. 4) 
c 

X 

P 

c 

n 

n 


<v< 

(c d } t 
c/R . 

t 


V 

8 


D 


(c a /c t> 

% 
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PROGRAM 

A7 

A8 

B 

B$ 

B(I) 

Bl 

B2 

B3 

B5 

B6 

B7 

C(1,I) 

C(2,I) 

C(3,I) 

C(4,I) 

C$ 

Cl 

C2 

C5 

C6 

C7 

C8 

C9 

D 

D(1,I) 

D(2,I) 

D(3,I) 

D(4,I) 

D$ 

D3 


USAGE 

Cumulative alpha change during off-design 
analysis, degrees 

Logical: Available thrust within 1% of 

required thrust (A8=0){ otherwise A8=l 

Number of propeller blades 

Output format string variable 

Beta at each blade station, degrees 

2 2 

Thrust coefficient, 2T/pV ifR 
Advance ratio 

Power Coefficient, 2P/PV 3 R 2 
Efficiency 

Blown propeller efficiency 

Logical: Beta fixed? l=yes, 0=no 

Lift Coefficient 

Drag Coefficient 

Chord/span ratio 

Unused 

Output format string variable 
Temporary variable 

Velocity of sound (f/s) at current altitude 
Case data: Chord fixed this run? l=yes, 0=no 

Case data: Momentum coefficients fixed this 

run? l=yes, 0=no 

Case data: Evaluate which propellers? 0=none, 

l=non-blown only, 2=blown only, 3= both 

Logical: Cross plot this run? l=yes, 0=no 

Temporary variable: maximum alpha, degrees 

Propeller diameter, feet 

Unused 

Unused 

Lookup table drag interval coefficients 
Drag to lift coef'icent ratio 
Output format string variable 
Temporary variable 




i 

i 






ENGINEERING 

PROGRAM 

USAGE 

^(Ref. 4) 

E(I) 

Vortex sheet spacing parameter 

^(Ref. 4) 

F(I) 

Ratio of average velocity increment in the 
slipstream to the sheet velocity 


FI 

Temporary variable 


FORl$ 

Output format string variable 


F0R2$ 

Output format string variable 


FOR3$ 

Output format string variable 


F0R4$ 

Output format string variable 


F0R5$ 

Output format string variable 


FTEM$ 

Output format string variable 

G . (Ref . 4) 

G(I) 

Circulation distribution function 

V 

(P T >. 

H(1,I) 

Jet velocity at each blade station, f/s 

H(2,I) 

Jet pressure at each blade station, PSF 


H(3,I) 

Local velocity at each blade station, f/s 

T 

H(4,I) 

Differential thrust; also total thrust, lbs. 

Q 

H(5,I) 

Differential torque; also total torque, ft-lbs. 

8 in 4. 
V 

H(6, I) 

Sin Phi-i 

cos 

H(7,I) 

Cos Phi-i 

• •> 
m 

H(8,I) 

Mass flow/foot, slugs/sec/foot 

HP ,• 
pet 

H(9,I) 

Horsepower/foot 

P. 

H 

Air density, slugs/f**3 

h 

H2 

Altitude, km 


I 

Program loop control variable; maximum 
value is number of blade stations 


18 

Program termination control variable 


J 

Program loop control variable 


J( ) 

Blown propeller lookup table interval 
lift and drag coefficients 


K 

Program loop control variable 

5 

K(I) 

Blade station radius to blade radius ratio 


K2 

Temporary variable 


LC1,I) 

Non-blown propeller tabl lookup interval 
lift coefficients 


L(2,I) 

Unused 


47 



i 

l 

I 


l 


f 


I 



l 

i 

| 


i 

i 





mir*w 





ENGINEERING 

PROGRAM 

USAGE 


M 

Temporary variable 

(HN) i 

M(I) 

Local blade station Mach number 

MO 

Total mass flow, slugs/second 


N 

Temporary variable 

n 

N1 

Engine speed, revolutions /second 


N2 

Engine speed change/ iteration required to 
match required thrust, revolutions/second 


N3 

Cumulative engine speed change in matching 



thrust, revolutions/second 


0 

Temporary variable 


P 

Temporary variable 

HP total 

PI 

HP; Engine power’ als0 » Required power 

P co 

P2 

Static Pressure, lbs/f**2 

n 

P3 

Pi (3.14159...) 


Q 

Temporary variable 

V co 

Q(Case#,0) 

Airspeed, f/s 

n 

Q(Case#,l) 

Engine speed, rev. /sec 

D 

Q(Case#,2) 

Propeller diameter, feet 

T 

Q(Case#,3) 

Thrust , lbs . 

h 

Q(Case#,4) 

Altitude, feet 

B 

Q(Case#,5) 

Number of propeller blades 

HP 

avail 

Q(Case#,6) 

Engine power, HP 


Q(Case#,7) 

Logical* Alpha fixed during thrust matching? 
l=yes, 0=no 


Q(Case#,8) 

Logical: Beta fixed during thrust matching? 

l=yes, 0=no 


Q(Case#,9) 

Logical: Momentum coefficients fixed during 

thrust matching? l=yes, 0=n®. 


Q(Case#,10) 

Propeller type selection: l=non-blown, 

2=blown, 3= both 

oil 

Q(Case#,ll) 

Alpha at blade station 1, degrees 

oia 

Q( Case#, 12) 

Alpha at blade station 2, degrees 

<*20 

• 

Q(Case#,30) 

• 

Alpha at blade station 20, degrees 


48 


— tfKMOte ' 






ENGINEERING PROGRAM 


Q 


Q6 

<“)< 


R(I) 



R 

R 


R6 

a . (Ref . 

4) 

T(1,I) 

(Ref. 

4) 

T(2,I) 

^ oi " 


T(3,I) 

a oi ” 


T(4,I) 

h 


T(5,I) 

h 


T(6,I) 

h 


T(7,I) 

h 


T(8,I) 

h 


T(9,I) 

h 


T(10,I) 

(C/R)^ 


T(11,I) 

(C/R) i 


T(12,I) 

(C/R) i 


T(13,I) 

(C/R)^ 


T(14,I) 

(C/R)^ 


T(15,I) 

(C/R)^ 


T(16,I) 

T 


T 

HP 


TO 

c 



HP 


T1 

pc 



T 


T2 



T3 

HP 


T4 

aero 


HP 


T5 

total 


T 


T6 



T7 

Case 


T8 

HP .. 
avail 


T9 

(C ) . 
'p 


U(1,I) 


USAGE 

Torque, ft-lbs 

Reynolds number at each blade station 
Temporary variable 
Propeller tip radius, feet 
Induced velocity equations parameter a-i 
Induced velocity equations parameter a’-i 
Induced velocity equations corrected Phi 
Induced velocity equations corrected Alpha 
Final Beta for Case 1, analytic propeller (A*) 
Final Beta for Case 2 propeller (B) 

Final Beta for Case 3, analytic propeller (C') 
Final Beta for Case 4 propeller (D) 

Final Beta for Case 1, blown propeller (A) 

Final Beta for Case 3, blown propeller (C) 

Final Chord ratio for propeller A* 

Final Chord ratio for propeller B 
Final Chord ratio for propeller C' 

Final Chord ratio for propeller D 
Final Chord ratio for propeller A 
Final Chord ratio for propeller C 
Thrust , lbs . 

Compressor power required to pump from P-static 

to required hub jet pressure, HP 

Total centrifugal horsepower requirement, HP 

Ambient temperature, Degrees R 

Total useful work, ft-lbs 

Aerodynamic torque, HP 

Total power required, HP 

Thrust developed, lbs. 

Logical t !=both blown & non-blown evaluation, 

0- either blown or non-blown evaluation 
Case identification* 0=manual input, 

1- 10 predefined in data statements 
Available horsepower @ current engine speed, HP 
Momentum Coeffi cents at each blade station 
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ENGINEERING 

PROGRAM 

USAGE 


U(2-7,I) 

Not used 

♦< 

a. 

i 

e 

U(8,I) 

Phi at each blade station, degrees 

U(9,I) 

Alpha at each blade station, degrees 

U(10,I) 

K(I) in the Induced velocity iterative equation 

a 

UO 

Angle of Attack in lift/drag lookup subroutine 


U1 

Temporary (lookup subroutine) 


U2 

Viscosity of air 


U4 

Temporary variable 

< P A 

U5 

Total hub pressure, lbs/ft**2 

h 

U6 

Temporary variable (Induced velocity equations) 


U9 

Temporary variable CATAN argument) 

v * 

V 

Free stream velocity, f/s 

V i 

V(I) 

Local velocity at each blade station, f/s 


VI 

Temporary variable (Induced velocity equations) 


V2 

Temporary variable (Induced velocity equations) 


V3 

Identification number of first blade station 
used in calculation of non-blown propeller 
lift and drag coefficients 


V4 

Logical: 0=Blown, l=Non- blown 


V7 

Logical: Induced velocity iterations complete? 

l=yes, 0=no. 


V8 

Logical: Diagonistic printer output desired? 

l=yes, 0=no 


V9 

Logical: Analytic output complete? l=yes.. 0=no 


W2 

Logical: Want CRT iterative output? l=yes 


X(I) 

x-i, 0mega*R/V 


Y 

Temporary variable 


Z(I) 

Simpson's rule transfer parameter array. (Value 
of integrand at each station.) 


Z 

Simpson's rule integral (Integral of Z(I)) 

C (Ref. 4) 

zo 

Displacement Velocity Ratio 


Z9 

Degrees/Radian conversion factor. 


w 


v 






50 



ORIGINAL PAGE 19 
OF POOR QUALITY 




Program Listing and Cross References 


10 

20 REM* * 

30 REM* PROPELLOR DESIGN PROGRAM FOR USE IN * 

40 REM* THE DESIGN OF A FULL-BLOWN PROP OR A NON-BLOWN PROP * 

50 REM* WHB LARC 1/18/83 * 

55 REM* * 

70 DIM E(20) ,X(20) ,F(20) ,G(20) ,R(20) ,H( 9,20) ,T( 16,20) ,Q( 10, 30) 

80 DIM M(20) ,K(20) ,U(10,20), J(9, 8) 

90 DIM L(2,20),D(4,20),A(4,20),Z(20),C(4,20),B(20),V(20) 

100 REM 

110 REM START INPUT OF VALUES 

120 

130 OPEN "OUTPUT.DAT" FOR OUTFUT AS FILE #1 % 

140 MARGIN #1%, 132% 

150 OPEN "CRPLOT.DAT" FOR OUTPUT AS FILE #2% 

160 MARGIN #2%, 80% 

170 DATA -.68, 53. 5,. 26, -1.0,. 009, -.40, -.001,0. 

180 DATA -.42,52.5, .27, -1.0, .008, -.40,0. ,0. 

190 DATA -.15, 51. 5,. 28, -2.0,. 008, -.40,. 0009, -.025 
200 DATA .13,49.5, .25, -.50,. 0089, -.425,. 0011,. 02 
210 DATA .38, 49.0,. 28, -3. 50,. 01, -.405,. 0031,. 04 
220 DATA .66, 45. 5,. 26, -5. 50,. 0131, -.36 c ,. 0024,. 075 
230 DATA . 92, 40.0,. 24, 08. 00,. 0155, -.29,. 0069,. 0 
240 DATA 1.16, 32.0,. 17, -9. 50,. 0224, -.29,. 0035,. 02 
250 DATA 1.33, 22. 5,. 07, -5. 50,. 0259, -.27,. 0077,. 27 
260 FOR 1-1 TO 9 
270 FOR J=1 TO 8 
280 READ J(I,J) 

290 NEXT J 
300 NEXT I 

310 DATA 270., 41. 6, 6., 324., 10000., 3., 0., 0,1, 0,3, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4 
320 DATA 4, 4, 4, 4, 4 

330 DATA 270. ,41.6,6,324. ,10000. ,3. ,0,1, 1,0, 2, -12. ,-11.49,-10. 92,-10.38,-9.82 
340 DATA -9. 29, -8. 7, -8. 12, -7 i, -6. 76, -5. 89, -4. 9,-3. 82,-2. 72 
350 DATA -1.55, -.3, 1.0, 2. 3, 3. 62, 5. 

"60 DATA 270., 41. 6, 6., 324., 10000., 3., 0., 0,1, 0,3, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2 
370 DATA 2, 2, 2, 2, 2 

380 DATA 270, 41. 6, 6, 324, 10000, 3, 0,1, 1,0, 2, -12., -11. 49, -10. 92, -10. 38, -9. 82 
390 DATA -9. 29, -8. 7, -8. 12, -7. 5, -6. 76, -5. 89, -4. 9, -3. 82, -2. 72 
400 DATA -1.55, -.3, 1.0, 2. 3, 3. 62, 5. 

410 DATA 125,45,6,347,0,3,0,1,1,1,2,4,4,4,4,4,4,4,4,4,4,4,4,4,^,4,4,4,4,4,4 
420 DATA 125, 45, 6, 347, 0,3, 0,1, 1,1, 2, -12., -11. 49, -10. 92, -10. 38, -9. 82 
430 DATA -9. 29, -8. 7, -8. 12, -7. 5, -6. 76, -5. 89, -4. 9, -3. 82, -2. 72 
440 DATA -1.55, -.3, 1.0, 2. 3, 3. 62, 5. 

450 DATA 125,45,6,347,0,3,0,1,1,1,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2 
460 DATA 125, 45, 6, 347, 0,3, 0,1, 1,1, 2, -12., -11. 49, -10. 92, -10. 38, -9. 82 
470 DATA -9. 29, -8. 7, -8. 12, -7. 5, -6. 76, -5. 89, -4. 9, -3. 82, -2. 72 
480 DATA -1.55, -.3, 1.0, 2. 3, 3. 62, 5. 

490 DATA 125,45,6,347,0,3,0,1,0,1,1,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4 
500 DATA 125,45,6,347,0,3,0,1,0,1,1,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2 
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ORIGINAL PASS; Vi 

510 FOR 1=1 TO 10 OF POOR QUALITY 

520 FOR J=0 TO 30 
530 READ Q(I,J) 

540 NEXT J 
550 NEXT I 

560 REM CLEAR THE 'ANALYTIC OUTPUT COMPLETE* FLAG 

570 V9=0 

580 REM RESET THE THRUST MATCHING PARAMETERS 

590 A3=0 
600 A7=0 
610 N3=0 


620 PRINT DATE s "}DATE$(0) 

630 REM ENTER THE CASE IDENTIFICATION FOR THE CURRENT RUN 

640 PRINT "INPUT 'CASE #• (1-10) FOR THIS RUN, OR '0* FOR MANUAL INPUT "} 
650 INPUT T8 

660 IF T8<=0 THEN GOTO 810 
670 IF T8>10 THEN T8=10 
680 V=Q(T8,0) 

690 N1=Q(T8,1) 

700 D=Q(T8,2) 

710 T=Q(T8,3) 

720 H2=Q(T8,4) 

730 B=Q(T8,5) 

740 P1=Q(T8,6) 

750 C5=Q(T8,7) 

760 B7=Q(T8,8) 

770 C6=Q(T8,9) 

780 C7=Q(T8,10) 

790 IF T8< >0 THEN GOTO 990 

800 REM INPUT AIRSPEED 

810 PRINT "ENTER AIRSPEED , V , IN FT/SEC"; 

820 INPUT V 

830 REM INPUT ROTATIONAL SPEED 

840 PRINT "ENTER ROTATIONAL SPEED , N1 , IN RPS"; 

850 INPUT N1 

860 REM INPUT PROP DIAMETER 

870 PRINT "ENTER PROP. DIAMETER , D , IN FEET"} 

880 INPUT D 

890 REM INPUT THRUST 

900 PRINT "ENTER THRUST , T , IN POUNDS"; 

910 INPUT T 

920 REM INPUT ALTITUDE (IN FEET) 

930 PRINT "ENTER ALTITUDE , H , IN FEET"; 

940 INPUT H2 

950 REM ** *' r t **** * *** , A * ^*^ *^ ,<> ^^ ** ^ **^ ^ ** < ^*^ **^**** * **** 

960 REM FIRST SET OF INPUT VALUES ARE COMPLETE 

970 REM**************************************** 11 *********** 1 ******* 

980 REM Z9 = DEGREES/RADIAN 

990 P3=3. 1415927 
1000 Z9-180.0/P3 

1010 REM H2 IS CONVERTED TO METRIC UNITS (KM) 

1020 H2»(3.048E-04)*H2 
1030 IF H&47 THEN 1070 
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1040 
1C , 
1060 
1070 
1080 
1090 
1100 
1110 
1120 
1130 
1140 
1.50 
1160 
1170 
1180 
1190 
1200 
1210 
1220 
1230 
1240 
1250 
1260 
1270 
1280 
1290 
1300 
1310 
1320 
1330 
1340 
1350 
1360 
1370 
1380 
1390 

1400 

1401 

1402 

1403 

1404 

1405 

1406 

1407 

1408 

1409 

1410 
1420 
1430 
1440 
1450 
1460 
1470 


PRINT "ALTITUDE ENTERED AT SETUP TIME IS TOO LARGE" 

GO TO 930 

Rfci: ATMOSPHEPIC DENSITY SUBROUTINE IS CALLED 

GOSUB C930 

IP T8< >0 THEN GOTO 1210 
REM INPUT SECOND SET OF VALUES 

. ». j-.i. 

Ji 

REM INPUT NUMBER OF BLADES ON PROP 

PRINT "ENTER NUMBER OF BLADES ON PROP , B"; 

INPUT B 

REM INPUT HORSEPOWER 

PRINT "ENTER ENGINE HORSEPOWER , P"; 

INPUT PI 

REM * * ** * **************************************** * ** * **** 

REM SECOND SET OF INPUT VALUES ARE COMPLETE 

REM R6 = TIP RADIUS OF PROP 

R6=D/2 

REM B1 = THRUST COEFFICIENT <EQN. 15? 

B1=8*T/ ( h*V**2*D**2*P3 ) 

REM B2 = ADVANCE ANGLE <EQN. 9? 

B2=V/(P3*N1*D) 

REM B3 = POWER COEFFICIENT 

B3=2*P1*550 / (H*V**3*R6**2*P3) 

FOR 1=1 TO 20 

REM K(I) = RATIO OF RADIUS DISTANCE FROM I TO TIP <EQN. 17?-— 

K(I)=0.05*I 
NEXT I 

REM BLOWN OR NON- BLOWN DESIGN 

IF T8< ?0 THEN GOTO 1420 

REM MANUAL INPUT FOR BLOWN/NON-BLOWN EVALUATION ************** 

PRINT" ENTER '1' FOR NON-BLOWN PROPELLER DESIGN, '2* FOR BLOWN DESIGN," 

PRINT"OR '3' FOR BOTH "j 

INPUT C7 

IF C7<1 THEN C7=l 
IF C7 >3 THEN C7=3 
PRINT 

PRINT" ENTER CODE FOR PARAMETER YOU WISH TO VARY IN MATCHING THRUSTS." 
PRINT" ENTER 'l* FOR CHORD SCALING, *2' FOR BETA SCALING, AND » 

PRINT"' 3' FOR ALPHA SCALING. "; 

INPUT Cl 

IF CKl OR Cl >3 THEN GOTO 1402 
IF C1=1 THEN C5=0 ELSE C5=l 
IF Cl=2 THEN B7=0 ELSE B7=l 
IF Cl=3 THEN C6=0 ELSE C6=l 

REM ****** NOW SET UP CONTROL PARAMETERS FOR SELECTED EVALUATION ******** 
IF C7=2 THEN V4=0 ELSE V4=l 
IF C7>=3 THEN T7=l ELSE T7=0 
IF V4=0 THEN V9=l 


REM 


WANT DIAGONISTIC OUTPUT??? 1=YES, 0=N0 INTO V8 
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1480 

1490 

1500 

1510 

1520 

1530 

1540 

1550 

1560 

1570 

1580 

1590 

1600 

1610 

1620 

1630 

1640 

1650 

1660 

1670 

1680 

1690 

1700 

1710 

1720 

1730 

1740 

1750 

1/60 

1770 

1780 

1790 

1800 

1810 

1820 

1830 

1840 

1850 

1860 

1870 

1880 

1890 

1900 

1910 

1920 

1930 

1940 

1950 

1960 

1970 

1980 

1990 

2000 

2010 

54 


PRINT" ENTER 1 TO INCLUDE DIAGONISTIC PRINTER OUTPUT, "j 
PRINT "OTHERWISE ENTER ’O’.. "} 

INPUT V8 

IF V8; 0 THEN V8=0 
IF V8>1 THEN V8=l 

PRINT "ENTER *1' FOR INFLOW ITERATION CRT OUTPUT, OTHERWISE ‘O' 
INPUT W2 

IF W2>1 THEN W2=l 
IF W2? 0 THEN W2=0 
REM 

IF T8c >0 THEN GOTO 1790 
IF(T7=1) AND (V4=0) THEN GOTO 1820 


ORIGINAL PAGE Kf 
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MANUAL INPUT OF ALPHA AT EACH BLADE STATION 

A(1,I) = BLADE ANGLE OF ATTACK 
A(2,I) = INTERVAL CONSTANTS 
A(3,I) = INTERVAL INCREMENTAL ALPHA 


REM 
REM 
REM 
REM 
REM 
REM 

i=i 

PRINT "ALPHA(";I;")"; 

INPUT A1 

PRINT "HOW MANY VALUES 0F";A1; 

INPUT K 

FOR J=I TO K+I-l 
A(l, J)=A1 
IF J=20 THEN 1820 
NEXT J 
I=J+1 
GOTO 1690 
FOR 1=1 TO 20 
A(1,I)=Q(T8,I+10) 

NEXT I 

rkin 

REM NOW HAVE 20 VALUES FOR A(1,I) 

REM A(1,I)=20 ALPHA VALUES 

REM ASSIGN LOOKUP TABLE INTERVALS 

i 

A(2,l)=-3 

A(2,2)=0 

A(2,3)=3 

A(2,4)=6 

A(2,5)=9 

A(2,6)=12 

A(2,7)=14 

A(2,8)*17 

Aft I r k** 

REM NOW ASSIGN LIFT AND DRAG VALUES TO LOOKUP TABLE COEFFICIENTS 
R EM** ** * 

REM 

L(l,l)=0.001 
L(l,2)=0.37 
L(l,3)=0. 73 
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2020 

2030 

2040 

2050 

2060 

2070 

2080 

2090 

2100 

2110 

2120 

2130 

2140 

2150 

2160 

2170 

2180 

2190 

2200 

2210 

2220 

2230 

2240 

2250 

2260 

2270 

2280 

2290 

2300 

2310 

2320 

2330 

2340 

2350 

2360 

2370 

2380 

2390 

2400 

2410 

2420 

2430 

2440 

2450 

2460 

2470 

2480 

2490 

2500 

2510 

2520 

2530 


L(l,4)=1.01 

L(l,5;=1.3 

L(l,6)=1.52 

L(l,7)=1.66 

L(l,8)=1.66 

L(l,9)=0.001 

D(3,l)=0.0123 

D(3,2)=0.0081 

D(3,3)=0.01 

D(3,4)=0.0127 

D(3,5)=0.O181 

D(3,6)=0.0242 

D(3,7)=0.0365 
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REM WE NOW HAVE INTERVAL LIMITS AND LINEAR COEFFICIENTS FOR THE LIFT * 

REM AND DRAG COMPUTATIONS NOW COMPUTE LIFT * 

REM' 

IF V4=0 THEN GOTO 2540 
FOR 1=1 TO 20 

T(4,I)=A(1,I) 

NEXT I 
V3=l 

GOSUB 8410 

IF V9=l THEN GOTO 2540 
REM * ** * ** ** * * 

REM OUTPUT VALUES 

REM^^ f ** , '* ,w **'’^’ f **'^ Wf ****^********* , ”' f ^***'”>*** 

PRINT n ,CHR$(12%) !FORM FEED 

PRINT #1%,TAB(40); "ANALYTIC RESULTS FOR CASE #";T8 


TO PRINTER 


#1 

#1 

#i 

#i 

#i 

#i 

#i 


#i 

#i 

#i 

#i 

#i 

#i 

n 

#i 


V";V; ! FREE STREAM VELOCITY 

N";N1; ! ENGINE SPEED, REV/SEC 
D";D; ! PROPELLER DIAMETER, FEET 
T";T; ! REQUIRED THRUST, LBS 
RHO" ; H ! AIR DENSITY, RHO 

B"; B; ‘NUMBER OF PROPELLER BLADES 

H,KM";H2; ! ALTITUDE, KM 

P";Pl5 ! AVAILABLE ENGINE HORSEPOWER 

V/ND";V/(N1*D) ! ADVANCE ANGLE, LAMBDA, DEGREES 

LAMBDA" ;B2; ! ADVANCE RATIO 
TC"; Blj ! THRUST COEFFICIENT 
PC"; B3 ! POWER COEFFICIENT 


PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT #1 
PRINT #1 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 

PRINT #1 , "KSI CL 

PRINT 
FOR 1=1 TO 20 

PRINT #1 ,K(I),C(1,I),D(4,I),A(1,I) 
NEXT I 


D/L 


ALPHA" 


REM 


FIRST OUTPUT COMPLETE, COMPUTE FOR SECOND OUTPUT 
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2540 FOR 1=1 TO 20 

2550 REM ** ***** * ****** <EQN. 8^ ****i eit *****4e*i c k**ie* i c&M*' ir:: * * a ***** 
2560 E( I )=0 . 5*B*SQR( B2**2+l )*( 1-K( I ) ) / B2 

2570 REM ************** <EQN 5> A A A A A m A Ai* A jl , * M > Ar , iif , A ,, A , *<r*<F<p* 

2580 X(I)=K(I)/B2 

2590 REM ******* *** **** <EQN 7 > &4r&rit k4 r i r t r k & a a a A*A , & , A , ‘A , A i A # *& & AA ft ArArA 

2600 Y=EXP(-E(I)) 

2610 F(I)=2*ATN(SQR(1-Y**2)/Y)/P3 
2620 REM *** ** ** * ** ** ** <EQN 6 ^ 

2630 G(I)=F(I)*X(I)**2/(X(I)**2+1) 

2640 REM ************** <EQN 20 > INTEGRAND ******************* * **** 
2650 M(I)=4*K(I)*G(I)*(1-D(4,I)/X(I)) 

2660 NEXT I 
2670 FOR 1=1 TO 20 
2680 Z(I)=M(I) 

2690 NEXT I 
2700 GOSUB 8830 

2710 REM ************** <EQN 20 > INTEGRAL 
2720 M=Z 

2730 FOR 1=1 TO 20 

2740 REM AAAfftfrArArtfr AArfrAVrA <EQN 21 > INTEGRAND ************************ 
2750 Z(I)=M(I)/(2*(X(I)**2+1)) 

2760 NEXT I 
2770 GOSUB 8830 

2780 REM <EQN 21 > INTEGRAL 

2790 N=Z 

2800 FOR 1=1 TO 20 

2810 REM ************** <EQN 16 > INTEGRAND, FIRST PART irk&'kic & ick' S ck'tck 
2820 Z(I)=4*K(I)*G(I)*(1+D(4,I)*X(I)) 

2830 NEXT I 
2840 GOSUB 8830 

2850 REM ^'rVn'fVf ft &■ AArAr ■ ft ArAr <EQN 16 > INTEGRAL, FIRST PART ************* 
2860 0=Z 

2870 FOR 1=1 TO 20 

2880 REM ************** < EQN 16> INTEGRAND, SECOND PART *********** 
2890 za)=4*K(I)*G(I)*(l+D(4,I)*X(I))*X(I)**2/(2*(X(I)**2+l)) 

2900 NEXT I 
2910 GOSUB 8830 

2920 REM ************** < EQN 16> INTEGRAL, SECOND PART 
2930 P=Z 

2940 IF T=0 THEN 3070 
2950 REM 

2960 REM THE REQUIRED THRUST WAS SPECIFIED 

2970 REM 
2980 REM 

2990 REM AfAf Af ArAr ArArArAnY Ar 1 ArtYA? < EQN 19> ************** ^ 

3000 ZO=M*( 1-SQR( 1- ( 4*B1*N/M**2 ) ) ) / ( 2*N) 

3010 REM ************** < EQN 29> NUMERATOR 
3020 B3=0*Z0+P*Z0**2 

3030 REM < EQN 29> *********************** * ** ** *** * ** 

3040 B5=B1/B3 
3050 GOTO 3160 


«»*«■ 





■ i *£; 


ORIGINAL PAGE 13 
OF. POOR QUALITY 


REM***********^ 

REM THE ENGINE HORSEPOWER WAS SPECIFIED 

REM******************************************* * *** * *** * ******* 

Z0=O*(SQR(l-t(4*B3*P)/O**2)-l)/(2*P) 

B1=M*Z0-N*Z0**2 

REM **fnV**VrtV****VfVr <EQN 29 > 

B5=B1/B3 

REM BEGIN SECOND SET OF OUTPUTS 

IF V9=l THEN 3210 
PRINT #1 

PRINT #1 Z";Z0; 

PRINT #1 ETA" ; B5 

REM VnWnSrfc ** * is isis icicin'c REQUIRED TORQUE FOR CURRENT PROPELLER **** 
Q6=V**3*B1*H*D**2/ ( 16*N1*B5 ) 

REM *** hWrf vL JU AJUJtiA A REQUIRED HORSEPOWER FOR CURRENT PROPELLER * 
P1=2*P3*N1*Q6/ 550 

REM ^*****^****** THRUST PRODUCED BY CURRENT PROPELLER ****** 

T=2*P3*N1*Q6*B5/V 

IF V9=l THEN GOTO 3370 

PRINT #1 

PRINT #1 ," TORQUE" ;Q6; 

PRINT #1 ," HP"; PI; 

PRINT #1 ," THRUST" ;T 

REM***********^ f ************************************ * * * * * **************** 
REM** ** 

REM** BYPASS ANALYTIC CHORD AND BETA CALCULATIONS ** 

REM** IF OFF DESIGN EVALUATION (T8>4) ** 

REM** ** 

REM** isisisisisisisisisicisisisisisisisisisisisisisisisisisisisisisisisisisisisisi 'sisirkisisisitisisisisisisisisiSititititisisieieit'kitit 

FOR 1=1 TO 20 

REM ******* ****** < EQN 24> Vf**-***Vf****VfVf*Vf*****^V**********V?***Vf*****Yf*** 


IF T8< 5 THEN C( 3 , I)=4*P3*B2*G( I )*Z0/ ( B*SQR(X( I)**2+1)*C( 1,1)) 

REM ************* < EQN 2 5 B > ******************************************vr* 
IF T8< 5 THEN B(I)=Z9*ATN(B2*(l+Z0/2)/K(I))+A(l,I) 

REM ********* SAVE CALCULATED ANALYTIC VALUES FOR SUBSEQUENT CASE USE *** 
IF (T8=l) AND (V4=l) THEN T(5,I)=B(I) 

IF (T8=2) THEN C(3,I)=T(15,I) 

IF (T8=3) AND (V4=l) THEN T(7,I)=B(I) 

IF (T8=4) THEN C(3,I)=T(16,I) 

IF (T8=5) THEN B(I)=T(9,I) 

IF (T8=5) THEN C(3,I)=T(15,I) 

IF (T8=6) THEN B(I)=T(6,I) 

IF (T8=6) THEN C(3,I)=T(12,I) 

IF (T8=7) THEN B(I)=T(10,I) 

IF (T8=7) THEN C(3,I)=T(16,I) 

IF (T8=8) THEN B(I)=T(8,I) 

IF (T8=8) THEN C(3,I)=T(14,I) 

IF (T8=9) THEN B(I)=T(5,I) 

IF (T8=9) THEN C(3,I)=T(11,I) 

IF (T8=10) THEN B(I)=T(7,I) 

IF (T8=10) THEN C(3,I)=T(13,I) 

NEXT I 
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3600 IF V9=l THEN GOTO 3670 
3610 PRINT #1 

3620 PRINT #1 ,"C/R .5(C/R)C0SB .5(C/R)SINB B» 

3630 PRINT #1 
3640 FOR 1=1 TO 20 

3650 PRINT #1 ,C(3,I) , 0 . 5*C(3,I)*COS(B(I)/Z9) ,0. 5*C(3,I)*SIN(B(I)/Z9) ,B(I) 
3660 NEXT I 
3670 GOSUB 8930 

3680 REM VELOCITY OF SOUND (IN FT/SEC) 

3690 C2=SQR( 2403 . 0*T2) 

3700 REM VISCOSITY (LB*SEC/FT**2) 

3710 U2=(340.8+0.548*(T2-453.0))*(10**(-9» 

3720 FOR 1=1 TO 20 

3730 REM ********** LOCAL VELOCITY, F/S ************** * ********** 

3740 V( I ) =SQR( V**2+( 2*P3*K( I )*R6*N1 ) **2 ) 

3750 REM MACH NUMBER 

3760 M(I)=V(I)/C2 

3770 REM REYNOLDS NUMBER 

3780 R( I )=H*V (I )*C( 3 , I )*D/ ( 2*U2 ) 

3790 NEXT I 

3800 IF V9=l THEN GOTO 3870 
3810 PRINT #1 

3820 PRINT #1 , "MACH NO REYNOLDS NO" 

3830 PRINT #1 

3840 FOR 1=1 TO 20 

3850 PRINT #1 ,M(I),R(I) 

3860 NEXT I 


3870 

3880 

3890 

3900 

3910 

3920 

3930 

3940 

3950 

3960 

3970 

3980 

3990 

4000 

4010 

4020 

4030 

4040 

4050 

4060 

4070 

4080 

4090 

4100 

4110 


REM SET ’ANALYTIC OUTPUT COMPLETE’ FLAG TO 1 

V9=l 

REM*******---**** 

REM* CALCULATE EACH BLADE STATIONS LOCAL VELOCITY [EQN. 7] 

REM******************* 


REM 

FOR 1=1 TO 20 

H( 3 , I ) = ( V**2+( P3*K( I ) *D*N1 )**2 )** . 5 
NEXT I 

REM** 

REM** ROUTE PROGRAM FLOW BASED ON BLOWN OR NON-BLOWN OPTION 
REM** (BLOWN, V4=0; NON-BLOWN, V4=l) 

REM** 


REM ****************************************************************** 
REM********** SET THE JET VELOCITY AT THE TIP TO . 95 MACH *** ******* 
REM********** IF BLOWN AND THE DESIGN POINT EVALUATION ************* 
REM********** OTHERWISE, SHUT BLOWING OFF 


REM 

IF (V4=0) AND (T8< 5) THEN H(1,19)=.95*C2 ELSE H(l,19)=0. 

REM 

REM********** THE REQUIRED JET PRESSURE AT STATION 19 IS [EQN 2] ****** 


REM 

H(2 , 19)=0 . 5*H*H(1 , 19)**2.+P2 


REM 




* *SK 


St** 
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REM ********** THE HUB PRESSURE TO GIVE THIS STATION 19 IS [EQN 6] drft&ftdr 
REM 

U4= ( (P3*N1*K( 1 9)*D)**2 ) / ( 3432*T2 ) 

U5=H(2 , 19)/EXP(U4) 

REM 

REM* CALCULATE JET PRESSURE AT EACH BLADE STATION [EQN. 6] 

REM 

FOR 1=1 TO 19 

U4= ( (P3*N1*K( I )*D)**2 ) / ( 3432*T2 ) 

H(2,I)=U5*EXP(U4) 

REM 

REM* CALCULATE JET VELOCITY AT EACH STATION ALSO [EQN. 2] 

REM 

A5=(H(2,I)-P2)*2/H 
IF A5<=0 THEN A5=0. 

H(l, I)=A5**.5 
REM 

REM* * 

REM* CALCULATE EACH BLADE STATIONS LOCAL VELOCITY [EQN. 7] * 

REM* * 

H( 3 , I)=(V**2+(P3*K(I)*D*N1)**2)** . 5 
V( I )=SQR(V**2+( 2*P3*K( I )*R6*N1 )**2 ) 

REM 

f «*“*'* • ******** ***************** * ** 

REM* ... AND FINALLY, CALCULATE THE MOMENTUM COEFFICIENT 
REM* AT EACH BLADE STATION [EQN. 1] 

REM 

REM 

U(l,I)=((H(l,I)/H(3,I))**2)/696.96 

REM******************************************************************** 

REM* END OF LOOP 

REM 

NEXT I 

REM 

REM * RECALCULATE [EQN 7] FOR CORRECT LOCAL VELOCITY DURING THE * 

REM * THRUST MATCHING ITERATIVE PORTION OF THIS PROGRAM. * 

FOR 1=1 TO 20 

H( 3 , I)=(V**2+( P3*K( I )*D*N1 )**2 )** . 5 
V( I)=SQR(V**2+( 2*P3*K( I )*R6*N1 )**2 ) 

NEXT I 

REM*********************************************************************** 

REM 

REM MAKE INFLOW CORRECTIONS FOR ALPHA, PHI, & RECALCULATE CL, CD, & CD/ CL 
REM 

REM 


59 




Si , «/*l ,i U-ItViit. 




4660 

4670 

4680 

4690 

4700 

4710 

4720 

4730 

4740 

4750 

4760 

4770 

4780 

4790 

4800 

4810 

4820 

4830 

4835 

4840 

4850 

4860 

4870 

4880 

4890 

4910 

4920 

4930 

4940 

4950 

4960 

4970 

4980 

4990 

5000 

5010 

5020 

5030 

5040 

5050 

5060 

5070 

5080 

5090 

5100 

5110 

5120 

5130 

5139 

5140 

5141 

5142 

5143 
5150 
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REM 

REM DEFINE PHI (I) AND 

REM 

B2=V/(P3*Nl*D) ! 

FOR 1=1 TO 20 

U(8,I)=ATN(B2/K(I)) ! 

A(1,I)=B(I)-U(8,I)*180. /P3 ! 
T(4,I)=A(1, I) ! 

NEXT I 

T(2,20)=0. ! 

T(3,20)=0. ! 

T(4,20)=U(8,20) ! 

T(l,20)=0. ! 

A(l,20)=0. ! 


OTHER QUANTITIES TO START ITERATION 
LAMBDA 

PHI(I) [EQN Al] 

CALCULATED ALPHA [EQN A2] 

ALPHA(I), CORRECTED 

a' AT EA’H STATION 
PHI, CORRECTED 
ALPHA, CORRECTED 
a AT EACH STATION 
ALPHA AT BLADE TIP 


REM 

REM INITIALIZE ITERATIVE EQUATIONS 

REM 


#.####### 


V7=l 

REM [EQN A3] 

FOR 1=1 TO 19 
H(6,I)=SIN(U(8,I)) 

H(7 , I )=C0S(U(8 , I ) ) 

V1=(((B/2)*(B2**2+1.)**.5)/B2)*(1.-K(I)) 

U6=EXP(-V1) 

U9=ATN(SQR( 1-U6**2 ) /U6 ) 

U(10,I)=B*C(3,I)/(8*P3*K(I))*l./(2/P3*U9) 

NEXT I 
V3=l 

IF V4=0 THEN GOSUB 9200 ELSE GOSUB 8410 

for5$=" ## m jwjw ###.## m.m #»#.#» 

REM * LIHIT AJb?HA TO 17 DEC (NON- BLOWN) OR 14 DEG (BLOWN) * 

FOR 1=1 TO 19 

IF V4=0 THEN C9*14.0 ELSE C9=17.0 
IF I >19 THEN C9=17,0 

IF (V4=0) AND (I<20) THEN Cl=1.4 ELSE Cl=1.66 
IF A(1,I)>C9 THEN C8=C1 ELSE C8=C(1,I) 

REM ** * irk ***** [EQN A4] & [EQN A5] 

V1=C8*H( 7 , I ) /H( 6 , I )**2 
V2=C8/H(7, I) 

T(1,I)=(U(10,I)*V1)/(1.-U(10,I)*V1) 

T(2,I)=(U(10,I)*V2)/(1.+U(10,I)*V2) 

REM *********** [EQN A6] 

T(3,I)=ATN((B2/K(I;*((1+T(1,I))/(1-T(2,I))))) 

REM *********** [EQN A7) 

T(4,I)=A(1,I)+28.65*(U(8,I)-T(3,I)) 

REM *********** IF REQUESTED, OUTPUT CRT DIAGONISTIC DATA **** 

IF W2<>1 THEN GOTO 5150 
IF I>1 THEN GOTO 5143 
PRINT 

PRINT ”V4 I AOA-I AOA-C PHI-I PHI-C C-MU CHC3 -/RAD" 

PRINT USING F0R5$,V4,I,A(l,I),T(4,I),U(8,I)*Z9,(B(I)-T(4, i )),U(l,I),C(3,I) 
NEXT I 
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5170 REM ***** DIAGONSTIC PRINTER OUTPUT ****** 

5180 IF V8< >1 THEN GOTO 5270 

5190 PRINT #1Z," I AOAI AOAC PHII PHIC A(I) A' (I)” 

5200 FOR I®13 TO 15 

5210 PRINT HZ USING '## ##>### ##.### ##.### ##.### ##.### ##.####',& 

I,A(1,I),T(4,I),U(8,I)*180/P3,T(3,I)*180/P3,T(1,I),T(2,I) 

5220 NEXT I 

5230 REM ***** END OF DIAGONISTIC OUTPUT AT THIS POINT ******* 

5240 REM 

5250 REM TEST ALPHA-ALPHAC FOR CONVERGENCE 

5260 REM 

5270 FOR 1=1 TO 19 

5280 REM ******** [EQN A9] ************************************************ 

5290 IF ABS(A(1,I)-T(4,I))>0.5 THEN V7=0 
5300 IF V8< >1 THEN GOTO 5330 

5310 IF ABS(A(1,I)-T(4,I))>.5 THEN PRINT ”I,A(1,I),T(4,I) = "jl, A(1,I),T(4,I) 
5320 REM ********* [EQN A8] UPDATE CURRENT PHI AT EACH BLADE STATION ******* 
5330 U(8,I)=(B(I)-T(4,I))*P3/180. 

5340 REM ********* UPDATE CURRENT ALPHA GUESS ******** * **** * **************** 
5350 A(1,I)=T(4,I) 

5360 NEXT I 
5370 V3=l 

5380 REM *********** GEf NEW Cl & Cd VALUES FOR BLOWN OR NON-BLOWN PROPELLER*** 
5390 IF V4=0 THEN GOSUB 9200 ELSE GOSUB 8410 
5400 IF V7=0 THEN GOTO 4830 E^SE GOTO 5430 
5410 REM 

5420 REM END OF ITERATIVE SECTION. 

5430 REM 

5440 REM CALCULATE DIFFERENTIAL THRUST AT EACH STATION, I, AND 

5450 REM DIFFERENTIAL TORQUE AT EACH STATION 

5460 REM 

5470 FOR 1=1 TO 20 
5480 H(6,I)=SIN(U(8,I)) 

5490 H(7,I)=COS(U(8,I)) 

5500 REM ** VfVrSr * * */nV VnV * [EQN 9] DIFFERENTIAL THRUST 

5510 H( 4 , I )=0 . 5*H*V**2*( ( 1+T( 1 , I ) ) /H( 6 , I ))**2*B*C( 3 , I )*D/ 2 

5520 H(4,I)=H(4,I)*(C(1,I)*H(7,I)-C(2,I)*H(6,I)) 

5530 REM *Vf**Y*-Vf*>V***Vf [EQN 10] DIFFERENTIAL TORQUE UrMt * * *,V * ** tV&tofc VnWc 

5540 H( 5 , 1 )=D*K( I ) /4*H*V**2*( ( 1+T( 1 , 1 ) ) /H( 6 , 1 ) )**2*B*C( 3 , 1 )*D/ 2 
5550 H(5,I)=H(5,I)*(C(1,I)*H(6,I)+C(2,I)*H(7,I)) 

5560 Z(I)=H(4, I) 

5570 NEXT I 
5580 REM 

5590 REM INTEGRATE DT/DR TO FIND TOTAL THRUST, T6 

5600 REM 

5610 GOSUB 8830 

5620 T6=Z*D*0 . 5 

5630 REM ********** TEST FOR CONVERGENCE OF ACTUAL THRUST TO REQUIRED THRUST ** 
5640 IF ABS(T-T6X T/100. THEN A8=0 ELSE A8=l 
5650 IF (T8>4) OR (A8=0) THEN GOTO 6520 
5655 PRINT 
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5660 
5670 
5680 
5690 
5700 
5710 
5720 
5730 
5/40 
5750 
5760 
5770 
5780 
5790 
5800 
5810 
5820 
5830 
5840 
5850 
5860 
5870 
5880 
5890 
5900 
5910 
5920 
5930 
5940 
5950 
5960 
5970 
5980 
5990 
6000 
6010 
6030 
6040 
6050 
6060 
6070 
6080 
6090 
6100 
6110 
6120 
6130 
6140 
6J 50 
6160 
6170 
6180 
6200 


REM ******* * i 
REM ** 

REM ** s 

REM ** 0 

REM ** C 

REM ** B 

REM ** c 

REM ** A 

REM ** 
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SCALE ALPHA, BETA, OR MOMENTUM COEFFICIENTS DEPENDING 
ON CASE BEING RUN. . . THE CONTROL VARIABLES ARE AS FOLLOW: 

C5=l (CASES 2,4-10, & BLOWN CASES) DO NOT RESCALE CHORDS 
B7=l (CASES 1-8) DO NOT RESCALE BETA 
C6=l (CASES 1N,3N, 5-10) DO NOT RESCALE ALPHA 
ALL QUANTITIES ARE RESCALED BY REQUIRED THRUST/ AVAIL. THRUST** 

** 


REM PRINT DIAGONISTICS HERE IF V8=l **************** 

IF V8<>1 THEN GOTO 5860 

PRINT #1%, "REQUIRED THRUST = ";T;" AVAILABLE THRUST = ";T6 
FOR 1=13 TO 15 

IF C5< >1 THEN PRINT' if 1%, "ALPHA VALUES ARE "jI,A(l,I) 

IF B7< >1 THEN PRINT #1%, "BETA VALUES ARE ";I,B(I) 

REM 
NEXT I 

REM ********* END OF DIAGONISTIC OUTPUT AT THIS POINT 
IF B7=l THEN 6030 

REM * RESCALE BETA TO MATCH AVAILABLE THRUST TO REQUIRED THRUST 

A2=.005*(T-T6) 

PRINT"REQUIRED THRUST=";T;" AVAILABLE THRUST=" ; T6 ; "AUTODELTA BETA="jA2 
PRINT"T0TAL CHANGE IN BETA SO FAR THIS RUN IS ";A3 
PRINT"ENTER ’O' IF OK, 'l* TO ENTER MANUAL CHANGE IN BETA," 

PRINT "OR '-l* TO FORCE END OF RUN..."} 

INPUT A4 

IF A4=-l THEN GOTO 6520 
IF A4=0 THEN GOTO 6000 

PRINT "ENTER NEW CHANGE IN BETA (DEGREES) "; 

INPUT A2 
A3=A3+A2 

PRINT "TOTAL CHANGE IN BETA NOW IS ";A3 
IF (C6=l) OR (V4=l AND T8>0) THEN GOTO 6200 

REM * RESCALE ALPHA TO MATCH AVAILABLE THRUST TO REQUIRED THRUST * 
A6= . 035*( T-T6 ) 

PRINT"REQUIRED THRUST=";T;" AVAILABLE THRUST=" ; T6 } "AUTODELTA ALPHA=";A6 
PRINT"TOTAL CHANGE IN ALPHA SO FAR THIS RUN IS ";A7 
PRINT"ENTER 'O' IF OK, '1' TO ENTER MANUAL CHANGE IN ALPHA," 

PRINT "OR '-1' TO FORCE END OF RUN..."; 

INPUT A4 

IF A4=-l THEN GOTO 6520 

IF A4»0 THEN GOTO 6170 

PRINT "ENTER NEW CHANGE IN ALPHA 

INPUT A6 

A7*A7+A6 

PRINT"TOTAL CHANGE IN 'ALPHA' IS NOW ";A7 
IF (T8<5) OR (T8?6) THEN GOTO 6390 
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6210 REM ****************Vf*vnV**v?*******VfVf*******Vf****foV****Vf**********Vc******* 
6220 REM * RESCALE ENGINE RPM TO MATCH AVAILABLE THRUST TO REQUIRED THRUST * 
6230 REM *********y~v*************************************^^^^ 

6240 N2=0 . 015*(T6-T) 

6250 PRINT"REQUIRED THRUST=";T;" AVAILABLE THRUST=";T6;" AUTODELTA SPEED=";N2 

6260 PRINT "TOTAL CHANGE IN SPEED SO FAR THIS RUN IS ";N3 

6270 PRINT "ENTER ’O’ IF OK, 'l* TO ENTER MANUAL CHANGE IN SPEED," 

6280 PRINT "OR ’-1' TO FORCE END OF RUN.... "; 

6290 INPUT A4 

6300 IF A4=-l THEN GOTO 6520 
6310 IF A4=0 THEN GOTO 6340 

6320 PRINT "ENTER NEW CHANGE IN SPEED (RPS) "; 

6330 INPUT N2 
6340 N3=N3+N2 

6350 PRINT "TOTAL CHANGE IN SPEED IS NOW "}N3 
6360 N1=N1+N2 
6370 GOSUB 9640 

6380 PRINT "AVAILABLE HORSEPOWER IS ";T9 
6390 FOR 1=1 TO 20 


6400 REM y^***************************y^**********************ttV************* 


6410 REM * RESCALE CHORDS TO OBTAIN REQUIRED THRUST * 

6422 IF B7<J>1 THEN B(I)=B(I)+A2 

6424 IF (C6< >1) AND (V4=0) THEN B(I)=B(I)+A6 

6430 IF (C5< >1) AND (V4=l) THEN C(3,i)=T/T6*C(3,I) 

6431 IF (C5< >1) AND (T8=0) THEN C(3,I)=T/T6*C(3,:0 

6440 NEXT I 


6450 REM **************************************************************** 


6460 REM ** 


** 


6470 REM ** FOR SCALED VALUES, PROGRAM FLOW GOES TO INFLOW ** 

6480 REM ** ITERATIONS ** 


6490 REM ** 


** 


6500 REM ****************************** 


*********************************** 


6501 IF W2< >1 THEN GOTO 4560 

6502 IF (C5c->1) AND (V4=l OR T8=0) THEN PRINT" CHORDS HAVE BEEN RESCALED" 

6510 GOTO 4560 

6520 REM 

6530 REM INTEGRATE DQ/DR TO FIND TOTAL TORQUE, Q6 

6540 REM 

6550 FOR 1=1 TO 20 
6560 Z(I)=H(5,I) 

6570 NEXT I 
6580 GOSUB 8830 


6550 Q6=Z*D*0 . 5 
6600 REM 
6610 REM 
6620 REM 
6630 T0=0. 

6640 T1=0. 

6650 REM 
6660 REM 
6670 REM 

6680 IF V4=l THEN GOTO 7030 


SET BLOWING HP VALUES TO ZERO FOR NON-BLOWING CASE 


SKIP BLOWING CALCULATIONS FOR NON-BLOWN CASE (V4=l) 
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6690 

6700 

6710 

6720 

6730 

6740 

6750 

6760 

6770 

6780 

6790 

6800 

6810 

6820 

6830 

6840 

6850 

6860 

6870 

688 

6890 

6900 

6910 

6920 

6930 

6940 

6950 

696C 

6970 

6980 

6990 

7000 

7010 

7020 

7030 

7040 

7050 

7060 

7070 

7080 

7090 

7100 

7110 

7120 

7130 

7140 

7150 

7160 

7170 

7180 

7190 

7200 

7210 

7220 

7230 
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CALCULATE THE MASS FLOW PER FOOT (H(8,I)) AND THE 
HORSEPOWER PER FOOT (H(9,I)) REQUIRED... 
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REM 
REM 
REM 
REM 

FOR 1=1 TO 20 
IF I >19 THEN GOTO 6820 
IF H(1,I)>0. THEN GOTO 6790 
H(8,I)=0 
GOTO 6810 

REM [EQN 11] 

H(8,I)=U(1,I)*C(3,I)*D/4*H*H(3,I)**2/H(1,I) 

H(9,I)=H(8,I)*6006.*T2*((H(2,I)/U5)**.286-1)/550. 

IF I >19 THEN H(9,I)=0 
IF I >19 THEN H(8,I)=0. 

Z(I)=H(9,I) 

NEXT I 
REM 

REM INTEGRATE HP REQUIRED/FOOT TO FIND TOTAL HP, T1 

REM 

GOSUB 8830 
T1=Z*B*D*.05 
REM 

REM INTEGRATE MASS FLOW/FOOT TO FIND TOTAL MASS FLOW 

REM 

FOR 1=1 TO 20 
Z(I)=H(8,1) 

NEXT I 
GOSUB 8830 
M0=Z*B*D*0 . 5 
REM ** 

REM * CALCULATE COMPRESSOR HP REQUIRED [EQN. 13] 

REM 

T0=M0*6006 ./ 550 . *T2*( (U5/P2 )**. 286-1) 

REM 

REM * CALCULATE TOTAL USEFUL WORK [EQN. 15] * 

REM 

T3=T6*V/550. 

REM * CONVERT AERODYNAMIC TORQUE TO HORSEPOWER [EQN. A10] * 

REM * AND THEN FIND THE TOTAL HORSEPOWER REQUIRED [EQN. 14] * 

REM 

T4=2.*P3*Nl*Q6/550. 

T5=T0+T1+T4 
REM ******** *** ** **tf 
REM * CALCULATE EFFICIENCY [EQN. 16] * 

REM 

B6=T3/T5 
REM 

REM * CALCULATE LOCAL MACH NUMBER [EQN. All], REYNOLD’S NUMBER * 

REM * [EQN. A12] , AND DRAG-TO-LIFT RATIO [EQN. A13] * 

REM 

FOR 1=1 TO 20 
Q(0,I)=H(3,I)/C2 
R( I )=H*C( 3,I)*D/2*H(3,I)/U2 
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7240 D(4,I)=C(2,I)/C(1,I) 

7250 NEXT I 

7260 REM Mirk**^****** DETERMINE AVAILABLE HP AT THE CURRENT RPM 
7270 GOSUB 9640 

7310 REM * THE NEXT SECTION SAVES DATA FOR POST-PROCESS CROSS PLOTTING * 
7330 IF (T8<5) OR (A4=-l) OR (A8=0) THEN GOTO 7580 

7340 FORI $=”00 ## ###.## ###.## ###.## ###.## #####. ###.## ##.# ##.# ## ” 
7350 FTEM $«"#**#.## ###.####" 

7360 F0R1$=F0R1$+FTEM$ 

7370 F0R2$=”01 «# #«*.«### ###.#### ###.#### ###.### #.### ####.##" 

7380 F0R3$="02 «# #.## ###.#### ###.#### #####.## ##.#» #.### ##.### » 
7390 FTEM $="#.### «#######. M 
7400 F0R3$=F0R3$+FTEM$ 

7410 F0R4$=”03 ## ####.### #.####### ######.#" 

7420 PRINT #2% USING F0R1$, T8, A3, A7,V,N1, 112/(3. 048E-4),T2,P2,D,B,T6,B6 
7430 PRINT #2% USING F0R2$,T8,T0,T1,T5,T9,U5/P2,T 
7440 FOR I =1 TO 20 

7450 PRINT m USING F0R3$,T8,K(I),C(1,I),C(2,I),1./D(4,I),A(1,I),& 
C(3,I),B(I),Q(0,I),R(I) 

7460 NEXT I 

7470 IF V4=l THEN GOTO 7510 
7480 FOR I = 1 TO 20 

7490 PRINT #2% USING F0R4$,T8,H(1,I),U(1,I),H(2,I) 

7500 NEXT I 

7510 PRINT #2%, ”99 99 9999.999 9.9999999 9999.999” 

7520 GOTO 5780 

7540 REM ** ** 

7550 REM ** FORMATTED LINE PRINTER OUTPUT FOR ALL CASES ** 

7560 REM ** ** 

7580 A4=0 

7590 PRINT #1, CHR$(12%) 5F0RM FEED 

7600 PRINT #1 , TAB ( 60% ) ; "RESULTS FOR CASE #";T8 

7610 IF V4=l THEN GOTO 7640 ! NON-BLOWN CASE 

7620 PRINT #1,TAB(35%); "PROPELLER CHARACTERISTICS" ; TAB( 105%) ; & 

"JET CHARACTERISTICS" 

7630 GOTO 7650 

7640 PRINT #1,TAB(35%); "PROPELLER CHARACTERISTICS" 

7650 PRINT #1 

7660 PRINT #1," KSI LIFT COEF DRAG COEF L/D RATIO ALPHA 
"CHORD/RAD TWIST MACH NO. REYNOLDS"; 

7670 IF V4=l THEN 7690 

7680 PRINT #1 , TAB (101%) ; "JET VEL. MOM. COEF JET PRES."i 
7690 PRINT #1 
7700 FOR 1=1 TO 20 

7710 PRINT #1 USING ' #.## fU.HHi ###.##### ###.####*,& 

K(I) f C(1,I),C(2,I),1/D(4,I),A(1,I); 

7720 PRINT #1 USING • #.####### ##.### LM ########.#•,& 

C(3,I),B(I),Q(0, I),R(I)> 

7730 IF V4=l THEN 7750 
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7740 

7750 

7760 

7770 

7780 

7790 

7800 

7810 

7820 

7830 

7840 

7850 

7860 

7870 

7880 

7890 

7900 

7910 

7920 

7930 

7940 

7950 

7960 

7970 

7980 

7990 

8000 

8010 

8020 

8030 

8040 

8050 

8060 

8070 

8080 

8090 

8100 

8110 

8120 

8130 

8140 

8150 

8160 

8170 

8180 

8190 

8200 

9210 

8220 

8230 

8240 

8250 


PRINT #1 USING ' ####.#### #.####### ####.####' ,6. 

H(1,I),U(1,I),H(2,I)} 

PRINT #1 
NEXT I 

%S A 0««»L «« ■ 

print #1 of poor QOA-rn 

PRINT #1, TAB(37) ; "OPERATING CONDITIONS” 

PRINT #1 

A$=’ VEL®### . # F/S ENG SPD= ##.## RPS * 

B$= • ALT®###### . ## FEET AIR DEN®# . ###### SL/ CO 'T * 

C$='TEMP =####.## DEG-R AMB PRES®####.# PSF* 

D$=A$+B$+C$ 

PRINT #1 

PRINT #1 USING D$,V,N1,H2/ (3.048E-4) ,H,T2 ,P2 
PRINT #1 

PRINT #1, TAB (40) ; "PROPELLER DATA" 

PRINT #1 

A$“' DIA®##.# FT BLDS®## TH =####.## LBS • 

B$®’EFF =#.#### LAM®##.### RQD HP®####.## AVL HP®####.## • 

C$=A$-i-B$ 

PRINT #1 USING C$,D,B,T6,B6,B2,T5,T9 
IF V4=l THEN GOTO 8040 
PRINT #1 
PRINT #1 

PRINT #1, TAB(42) } "COMPRESSOR" 

PRINT #1 

A$=' MASS FLOW®###.#### SL/SEC COMP HP®###.#### • 

B$= ’ COMPRESSOR RATIO®### . #### ’ 

C$=A$+B$ 

PRINT #1 USING C$,M0,T0,U5/P2 
PRINT #1 


REM 
REM ** 
REM ** 
REM ** 
ri^M ** 


END OF OUTPUT FOR THE PRESENT CASE 


■ Me kltl r kl r ! ; ‘rki r kitit 

REM********************************* ******* * * * * * * ****************** ******* 
REM ** SAVE BETA, B(I), AND FINAL CHORDS, C(3,I) FOR PROPELLERS ** 
REM** A, A', B, C, C’, AND D ACCORDING TO SCHEME IN NOTES... ** 
REM*’ 


IF T8>4 THEN GOTO 8280 
FOR I®1 TO 20 

If (T8=l) AND (V4«l) THEN T(11,I)»C(3,I) 
IF (T8«l) AND (V4=0) THEN T(15,I)»C(3,I) 
IF (T8«l) AND (V4-0) THEN T(9,I)«BvI) 

IF (T8«2) THEN T(12,I)«C(3,I) 

IF (T8«2) THEN T(6,I)-B(I) 

IF (T8®3) AND (V4®1) THEN T(13,I)-C(3,I) 
IF (T8®3) AND (V4-0) THEN T(16,I)-C(3,I) 
IF (T8-3J AND (V4«0) THEN T(10,I)-B(I) 

IF (T8«4) THEN T(14,I)-C(3,I) 


it; 

*nt -J-< ■ 




Ar. 
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82o0 IF (T8=4) THEN T(8,I)=B(I) 

8270 NEXT I 

8280 IF (T7=0) OR (V4=0) THEN GOTO 8320 
8290 V4=0 
8300 GOTO 1790 

8310 REM kirkiciciricid kirk "EXAMINE NEXT CASE?" DECISION POINT •kkickkieici e irkir kirki rk i rk 

8320 PRINT "ENTER OPTION NUMBER" ;TAB( 30) ;"0 - STOP" 

8330 PRINT TAB(30);"1 = RESTART PROGRAM " 

8340 PRINT TAB ( 30 ) ; " CONTROL-C OR CONTROL-Y = EXIT TO MONITOR" 

8350 PRINT "ENTER OPTION..."; 

8360 INPUT 18 

8370 IF I8<=0 THEN STOP 

8380 IF I8>3 THEN STOP 

8390 GOTO 570 

8400 STOP 

8410 REM irkirkicirkiei ti r ki ei cifisic irk is i fisig irkirkirir!cirkiririrkirkirkirkiri% i r ki (icic , k i c , k'kiciciei:irkicicicicicicic 


8420 REM ** 


** 

8430 REM ** 

CALCULATE LIFT, DRAG, AND DRAG/LIFT RATIO FOR THE NON- 

** 

8440 REM ** 

BLOWN PROPELLER 

** 

8450 REM ** 

INPUTS: V3=STARTING LOOP VALUE (IxRST BLADE STATION) 

irk 

8460 REM ** 

T(4,I)=ALPHA*S 

irk 

8470 REM ** 

A( 2 , I )=INTERVAL CONSTANTS 

** 

8480 REM ** 

L(1,I)=LIFT CONSTANTS 

** 

8490 REM ** 

D(3,I)=DRAG CONSTANTS 

irk 

8500 REM ** 

OUTPUTS :C(1, I )=LIFT COEFFICIENTS 

irk 

8510 REM ** 

C(2,I)=DRAG COEFFICIENTS 

irk 

8520 REM ** 

D(4,I)=DRAG/LIFT RATIO 

irk 

8530 REM ** 


irk 


8540 REM **^*^****:WWHif^*^ 

8550 FOR I=V3 TO 20 
8560 D(3,8)=SIN(T(4, I)) 

8570 D(3, 9)= ;; IN(T(4,I)) 

8580 IF T(4,l. :=12 THEN Fl=6 ELSE Fl=7 
8590 IF T(4,I)>17. THEN GOTO 8720 
8600 REM 

8610 REM J= RANK INDEX FOR LIFT AND DRAG 

8620 REM 

8630 J=INT((T(4,I)+Fl)/3) 

8640 REM 

8650 REM A( 3, I) “DELTA ALPHA 

8660 REM 

8670 K2=L(1, J+l) 

8680 D3=A(2, J+l)-A(2, J) 

8690 C(l f I)=L(1, J)+( (K2-L(l , J))/D3)*(T(4,I)-A(2, J)) 

8700 C(2,I)=D( 5, J)+((D(3, J+l)-D(3, J))/D3)*(T(4,I)-A(2, J)) 

8710 GOTO 8740 

3720 0(1, I)-. 001 

8) 0 C(2,I)=SIN(T(4,I)/Z9) 

8740 IF ABS(C(2,I))< 0.0001 THEN C(2 ,I)=0.0001 
8750 D(4, I)=C(2,I)/C(1, I) 

8760 NEXT I 
8770 RETURN 



8 

I 

I 









A- ' 


_ V.I 


8780 

8790 

8800 

8810 

8820 

8830 

8840 

8850 

8860 

8870 

8880 

8890 

8900 

8910 

8920 


ftft* Aft ***** *** * ********************* 

REM ** ** 
REM ** SIMPSON’S RULE INTEGRATION ** 
REM ** ** 


Q=0 

FOR J=1 TO 10 
Q=Q+Z(1+2*(J-1)) 
NEXT J 
R=0 

FOR J=1 TO 10 
R=R+Z(2*J) 

NEXT J 

Z=0.05*(4*Q+2*R)/3 

RETURN 
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8930 REM 

8940 REM ** ** 

8950 REM ** ATMOSPHERIC CHARACTERISTICS SUBROUTINE ** 

8960 REM ** ** 

8970 REM Ic'kitic * £ £ ft ft ft ftft ft ft ftft ft ftft ft ft ft ft ftft ft ft ft ft ft ft ftftftft ftft ft ft ft ft -ft ft ft - ftft ftft ft ft 

8980 IF H2S41 THEN 9070 

8990 IF H2< 20 THEN 9100 

9000 IF H2< 32 THEN 9130 

9010 IF H2< 47 THEN 9160 

9030 REM H = DENSITY IN SLUGS/FT**3 

9040 REM P2 = AMBIENT PRESSURE IN LB/FT**2 

9050 REM T2 = AMBIENT TEMPERATURE IN DEGREES ’R’ 

9060 EEM **^^** *** * * ^*^^^^**^****** *** ^***^'^*******^*^*^* 
9070 P2=2116.6792*(288. 15/(288. 15-6. 5*H2))**(-5. 255876) 

9080 T2=(288 . 15-6 . 5*H2)*1 . 8 
9090 GOTO 9180 

9100 P2=472 . 78248*EXP(-0. 157688*(H2-1I) ) 

9110 T2=389.97 
9120 GOTO 9180 

9130 P2=114 . 37003*(216 . 65/ ( 216 . 65+(H2-20) ) )**34 . 163195 
9140 T2=(216.65+(H2-20))*1.8 
9150 GOTO 9180 

9160 P2=18 . 132812*(228 . 65/ ( 228 . 65+2 . 8*(H2-32) ) >**12 . 201141 
9170 T2=(288.65+2.8*(H2-32))*1.8 
9180 H=0 . 0005827*P2/T2 
9190 RETURN 


9200 REM ftftftftftftftftftft ft ft ftftft ftftftftftftft ftft ftftftftftft ftft ftftftftft ft ftftftftftftftftftftftft ft ftftftftftft 


9210 REM ** ** 
9220 REM ** BLOWN PROPELLER LINEARIZED LOOKUP SUBROUTINE ** 
9230 REM ** **** INPUTS**** ** 
9240 REM ** MOMENTUM COEFFICIENTS AT STATION I t 17(1, I) ** 
9250 REM ** BLADE ANGLE OF ATTACK AT STATION I, T(4,I) ** 
92 SO REM ** ** 


9270 REM ** ****************** * ****^^ ft ftftftftft ft ftft ft ftftftftftrft ftft 
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9280 REM ** LIFT COEFFICIENTS AT STATION I, C(1,I) ** 
9290 REM ** DRAG COEFFICIENTS AT STATION I, C(2,I) ** 
9300 REM ** DRAG/LIFT COEFFICIENTS AT STATION I, D(4,I) ** 
9310 REM ** ** 

9330 FOR 1=1 TO 19 
9340 GOSUB 9470 


9350 IF (T(4,I)c -12.0)AND (W2=l) THEN PRINT "I ,A0A,U0,U1="; I ,T(4,I) ,U0,U1 
9360 IF U0i>=14.0 THEN GOTO 9400 

9370 C(1,I)=J(U0,1)+J(U0,2)*U(1,I)+U1*(J(U0,3)+J(U0,4)*U(1,I)) 

9380 C(2,I)=J(U0,5)+J(U0,6)*U(1,I)+U1*(J(U0,7)+J(U0,8)*U(1,I)) 

9390 GOTO 9420 
9400 C(1,I)=0 

9410 C(2,I)=SIN(T(4,I)/Z9) 

9420 IF ABS(C(1,I)X 0.001 THEN C(1,I)=0.001 
9430 IF ABS(C(2,I))c 0.0001 THEN C(2,I)=0.0001 
9440 D(4,I)=C(2,I)/C(1,I) 

9450 NEXT I 

9460 RETURN 


9480 REM** ** 

9490 REM** SUBROUTINE TO CALCULATE LOOKUP TABLE INTERVAL ** 
9500 REM** RANK UO AND AOA FRACTION U1 ** 

9510 REM** ** 

9520 REM* 

9530 U0=T(4,I)-.0001 
9540 IF UO -12. THEN 9580 
9550 IF U0>14.0 THEN 9620 
9560 U0=INT((U0+12)/3.0)+1.0 
9570 GOTO 9590 
9580 U0=1 

9590 Ul=(T(4,I)-(3*U0-15))/3.0 



9600 IF U0=9 THEN U1=U1*1.5 
9610 GOTO 9630 
9620 U'1~0 
9630 RETUx<N 

9640 REM ****’V******^Wf*5V**5V**5V»v*****Vf***Vf******?v****************** 


9650 REM ** ** 
9660 REM ** SUBROUTINE TO CALCULATE AVAILABLE HORSEPOWER ** 
9670 REM ** AS A FUNCTION OF ENGINE RPM @ FULL THROTTLE ** 
9680 REM ** ** 


9690 REM **’^*************>ViVVnV*VfVnV***********V.-Vf*****************i(lrthV 

9700 T9=( - . 1667)*Nl**2+20 . 8332*N1-314 . 99 
9710 RET'JRN 
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Symbol Cross Reference Table 

The following table lists for each program symbol 'variable name) the 
number (label) of every statement in which that symbol can be found. Symbols 
that end with the dollar sign ($) are string variables, and are used here to 
store formatting information for printer, terminal, or disk file output. All 
other symbols are used for numDeric data, with no differentiation made for 
integer or floating-point data. The Class column contains only the letter 'I* 
or a blank; if an 'I' is present, the corresponding Symbol is an array 
variable, otherwise, it is a scalar variable. 

The 'References' column lists the statement labels for those statements in 
which each symbol is found. The statement label is the integer portion of the 
Reference number; the fractional portion indicates which element of the 
statement contains the symbol, and is most often '.001' indicating the first 
part. The fractional portion will be other values when the symbol is found in 
the 'IF. . .THEN. . .ELSE. .. ' type statement, with the occurrance found aft«r 
'THEN' or 'ELSE'. The pound sign (#) indicates the statement in which array 
variables are dimensioned. 


Symbol 

Class 


References 



A 

I 90.001# 

1740.001 

1800.001 

1870.001 

1880.001 


1890.001 

1900.001 

1910.001 

1920.001 

1930.001 


1940.001 

2210.001 

2490.001 

3410.002 

4720.001 


4730.001 

4790.001 

5030.001 

5120.001 

5143.001 


5210.001 

5290.001 

5310.001 

5310.002 

5350.001 


5810.002 

8700.001 

7450.001 

7710.001 

8680.001 

8690.001 

A$ 

7820.001 

8020.001 

7850.001 

7910.001 

7930.001 

8000.001 

A1 

1700.001 

1710.001 

1740.001 



A2 

5900.001 

5910.001 

5990.001 

6000.001 

6422.002 

A3 

590.001 

5920.001 

6000.001 

6010.001 

7420.001 

A4 

5950.001 

5960.001 

5970.001 

6120.001 

6130.001 


6140. r, 01 
7580.001 

6290.001 

6300.001 

6310.001 

7330.001 

A5 

4290.001 

4300.001 

4300.002 

4310.001 


A6 

6070.001 

6080.001 

6160.001 

6170.001 

6424.002 

A7 

600.001 

6090.001 

61V j.001 

6180.001 

7420.001 

A8 

5640.002 

5640.003 

5650.001 

7330.001 
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Symbol 

Class 



References 



B 


730.001 

1140.001 

2370.001 

2560.001 

3390.002 



4870.001 

4910.001 

5510.001 

5540.001 

6900.001 



6980.001 

7420.001 

7940.001 



B 

I 

90.001# 

3410.002 

3430.002 

3450.002 

3470.002 



3490.002 

3510.002 

3530.002 

3550.002 

3570.002 



3650.001 

4720.001 

5143.001 

5330.001 

5820.002 



6422.002 

6424.002 

7450.001 

7720.001 

8190.002 



8210.002 

8240.002 

8260.002 



B$ 


7830.001 

8020.001 

7850.001 

7920.001 

7930.001 

8010.001 

B1 


1240.001 

2430.001 

3000.001 

3040.001 

3100.001 



3120.001 

3210.001 




B2 


1260.001 

2420.001 

2560.001 

2580.001 

3390. 0J2 



3410.002 

7940.001 

4690.001 

4710.001 

4870.001 

5100.001 

B3 


1280.001 

3120.001 

2440.001 

3020.001 

3040.001 

3090.001 

B5 


3040.001 

3120.001 

3190.001 

3210.001 

3250.001 

B6 


7160.001 

7420.001 

7940.001 



B7 


760.001 

6422.001 

1408.002 

1408.003 

5820.001 

5860.001 

C 

I 

90.001# 

2490.001 

3390.002 

3440.002 

3460.002 



3480.002 

3500.002 

3520.002 

3540.002 

3560.002 



3580.002 

3650.001 

3780.001 

4910.001 

5030. 0C3 



5143.001 

5510.001 

5520.001 

5540.001 

5550. C01 



6430.002 

6431.002 

6790.001 

7230.001 

7240.001 



7450.001 

7710.001 

7720.001 

8170.002 

8180.002 



8200.002 

8220.002 

8230.002 

8250.002 

8690.001 



8700.001 

8720.001 

8730.001 

8740.001 

8740.002 



8750.001 

9370.001 

9380.001 

9400.001 

9410.001 



9420.001 

9420.002 

9430.001 

9430.002 

9440.001 

C$ 


7840.001 

8030.001 

7850.001 

7930.001 

7940.001 

8020.001 

Cl 


1405.001 

1406.001 

1407.001 

1408.001 

1409.001 



5020.002 

5020.003 

5030.002 



C2 


3690.001 

3760.001 

4060.002 

7220.001 


C5 


750.001 

1407.002 

1407.003 

5810.001 

6430.001 



6431.001 

6502.001 




C6 


770.001 

1409.002 

1409.003 

6030.001 

6424.001 
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C7 

780.001 

1400.002 

1380.001 

1420.001 

C8 

5030.002 

5030.003 

C9 

5000.002 

5000.003 

D 

700.001 

2330.001 

4140.001 

5510.001 

6900.001 

880.001 

2400.001 

4220.001 

5540.001 

6980.001 

D 

I 90.001# 

2120.001 
2820.001 
8560.001 

2080.001 

2130.001 

2890.001 

8570.001 

D$ 

7850.001 

7870.001 

D3 

8680.001 

8690.001 

E 

I 70.001# 

2560.001 

F 

I 70.001# 

2610.001 

FI 

8580.002 

8580.003 

F0R1$ 

7340.001 

7360.001 

F0R2$ 

7370.001 

7430.001 

F0R3$ 

7380.001 

7400.001 

F0R4$ 

7410.001 

7490.001 

F0R5$ 

4950.001 

5143.001 
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1390.001 

1430.001 

1390.002 

1400.001 

5050.001 

5060.001 


5010.002 

5030.001 


1220.001 

3210.001 

4380.001 

5620.001 

7230.001 

1240.001 

3780.001 

4570.001 

6590.001 

7420.001 

1260.001 

3940.001 

4690.001 

6790.001 

7940.001 

2090.001 

2140.001 

7240.001 

8700.001 

2100.001 

2490.001 

7450.001 

8750.001 

2110.001 

2650.001 

7710.001 

9440.001 


8700.001 

2600.001 

2630.001 

8630.001 

7420.001 

7450.001 


7350.001 

7360.001 

7390.001 

7400.001 


70.001# 

3390.002 

2630.001 

2650.001 

2820.001 

2890.001 

70.001# 

4150.001 

4470.001 

5060.001 

5540.001 

6760.001 

6840.001 

7740.001 

3940.001 

4230.001 

4570.001 

5480.001 

5550.001 

6790.001 

6950.001 

4060.002 

4290.001 

4850.001 

5490.001 

5560.001 

6810.001 
7220.001 

4060.003 

4310.001 

4860.001 

5510.001 

6560.001 

6820.002 
7230.001 

4100.001 

4380.001 

5050.001 

5520.001 

6750.001 

6830.002 
7490.001 
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H 

1240.001 

1280.001 


4100. 001 

4290.001 


7230.001 

7870.001 

H2 

720.001 

940.001 


7420.001 

7870.001 


9010.001 

9070.001 


5140.001 

9160.001 

I 

260.001 

280.001 


550.001 

1290.001 


1690.001 

1730.001 


1810.001 

2200.001 


2490.001 

2500.001 


2600.001 

2610.001 


2670.001 

2680.001 


2760.001 

2800.001 


2890.001 

2900.001 


3430.002 

3440.002 


3480.002 

3490.002 


3530.002 

3540.002 


3580.002 

3590.001 


3720.001 

3740.001 


3840.001 

3850.001 


3950.001 

4210.001 


4310.001 

4380.001 


4560.001 

4570.001 


4710.001 

4720.001 


4850.001 

4860.001 


4990.001 

5010.001 


5050.001 

5060.001 


5120.001 

5140.001 


5210.001 

5220.001 


5310.002 

5330.001 


5480.001 

5490.001 


5550.001 

5560.001 


5820.002 

5840.001 


6430.002 

6431.002 


6570.001 

6730.001 


6790.001 

6810.001 


6830.002 

6840.001 


6960.001 

7210.001 


7250.001 

7440.001 


7490.001 

7500.001 


7740.001 

7760.001 


8190.002 

8200.002 


8240.002 

8250.002 


8560.001 

8570.001 


8690.001 

8700.001 


8740.002 

8750.001 


9350.002 

9370.001 


9420.001 

9420.002 


9450.001 

9530.001 
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QUALITY 

2350.001 

3210.001 

3780.001 

5510.001 

5540.001 

6790.001 

9180.001 

1020.001 

1030.001 

2380.001 

8980.001 

8990.001 

9000.001 

9080,001 

9100.001 

9130.001 

9170.001 

300.001 

510.001 

530.001 

1310.001 

1320.001 

1680.001 

1770.001 

1790.001 

1800.001 

2210.001 

2220.001 

2480.001 

2540.001 

2560.001 

2580.001 

2630.001 

2650.001 

2660.001 

2690.001 

2730.001 

2750.001 

2820.001 

2830.001 

2870.001 

3370.001 

3390.002 

3410.002 

3450.002 

3460.002 

3470.002 

3500.002 

3510.002 

3520.002 

3550.002 

3560.002 

3570.002 

3640.001 

3650.001 

3660.001 

3760.001 

3780.001 

3790.001 

3860.001 

3930.001 

3940.001 

4220.001 

4230.001 

4290.001 

4390.001 

4470.001 

4510.001 

4580.001 

4590.001 

4700.001 

4730.001 

4740.001 

4840.001 

4870.001 

4910.001 

4920.001 

5020.001 

5030.001 

5030.003 

5070.001 . 

5080.001 

5100.001 

5143.001 

5150.001 

5200.001 

5270.001 

5290.001 

5310.001 

5350.001 

5360.001 

5470.001 

5510.001 

5520.001 

5540.001 

5570.001 

5800.001 

5810.002 

6390.001 

6 422.002 

6424.002 

6440.001 

6550.001 

6560.001 

6740.001 

6750.001 

6760.001 

6820.001 

6820.002 

6830.001 

6850.001 

6940.001 

6950.001 

7220.001 

7230.001 

7240.001 

7450.001 

7460.001 

7480.001 

7700.001 

7710.001 

7720.001 

8160.001 

8170.002 

8180.002 

8210.002 

8220.002 

8230.002 

8260.002 

8270.001 

8550.001 

8580.001 

8590.001 

8630.001 

8720.001 

8730.001 

8740.001 

8760.001 

9330.001 

9350.001 

9380.001 

9400.001 

9410.001 

9430.001 

9430.002 

9440.001 

9590.001 
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8360.001 

8370.001 

J 


270.001 

280.001 



540.001 

1730.001 



1770.001 

8630.001 



8700.001 

8840.001 



8890.001 

8900.001 

J 

I 

80.001# 

280.001 

K 


1720.001 

1730.001 

K 

I 

80.001# 

1310.001 



2650.001 

2820.001 



3940.001 

4140.001 



4570.001 

4580.001 



5100.001 

5540.001 

K2 


8670.001 

8690.001 

L 

I 

90.001# 

1990.001 



2030.001 

2040.001 



8670.001 

8690.001 

M 


2720.001 

3000.001 

M 

I 

80.001# 

3850.001 

2650.001 

MO 


6980.001 

7020.001 

N 


2790.001 

3000.001 

N1 


690.001 

850.001 



3210.001 

3230.001 



4140.001 

4220.001 



4580.001 

4690.001 



7870.001 

9700.001 

N2 


6240.001 

6250.001 

N3 


610.001 

6260.001 

0 


2360.001 

3020.001 

P 


2930.001 

3020.001 

PI 


740.001 

3290.001 

1170.001 

P2 


4100.001 

4290.001 



7870.001 

8030.001 



9160.001 

9180.001 
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References 

8380.001 

cmi i 
QE POOR 

QUALITY 

290.001 

1740.001 

8670.001 

8850.001 

520.001 

1750.001 

8680.001 
8860.001 

530.001 

1760.001 

8690.001 

8880.001 

9370.001 

9380. 0C1 


2490.001 

2890.001 

4220.001 

4710.001 

7450.001 

2560.001 

3410.002 

4380.001 

4870.001 

7710.001 

2580.001 

3740.001 

4390.001 

4910.001 

2000.001 

2050.001 

2010.001 

2060.001 

2020.001 

2070.001 

3100.001 



2680.001 

2750.001 

3760.001 

8030.001 



3100.001 



1260.001 

3250.001 

4380.001 

6360.001 

2320.001 

3740.001 

4390.001 

7110.001 

2400.001 

3940.001 

4570.001 

7420.001 

6330.001 

6340.001 

6360.001 

6340.001 

6350.001 


3090.001 



3090.001 



1280.001 

2390.001 

3230.001 

7020.001 

9070.001 

7420.001 

9100.001 

7430.001 

9130.001 
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P3 


990.001 

1000.001 

1240.001 

1260.001 

1280.001 



2610.001 

3230.001 

3250.001 

3390.002 

3740.001 



3940.001 

4140.001 

4220.001 

4380.001 

4390.001 



4570.001 

4580.001 

4690.001 

4720.001 

4910.001 



5210.001 

5330.001 

7110.001 



Q 


8830.001 

8850.001 

8910.001 



Q 

I 

70.001# 

530.001 

680.001 

690.001 

700.001 



710.001 

720.001 

730.001 

740.001 

750.001 



760.001 

770.001 

780.001 

1800.001 

7220.001 



7450.001 

7720.001 




Q6 


3210.001 

3230.001 

3250.001 

3280.001 

6590.001 



7110.001 





R 


8870.001 

8890.001 

8910.001 



R 

I 

70.001# 

3780.001 

3850.001 

7230.001 

7450.001 



7720.001 





R6 


1220.001 

1280.001 

3740.001 

4390.001 

4580.001 

T 

I 

70.001# 

2210.001 

3430.002 

3440.002 

3450.002 



3460.002 

3470.002 

3480.002 

3490.002 

3500.002 



3510.002 

3520.002 

3530.002 

3540.002 

3550.002 



3560.002 

3570.002 

3580.002 

4730.001 

4750.001 



4760.001 

4770.001 

4780.001 

5070.001 

5080.001 



5100.001 

5120.001 

5143.001 

5210.001 

5290.001 



5310.001 

5310.002 

5330.001 

5350.001 

5510.001 



5540.001 

8170.002 

8180.002 

8190.002 

8200.002 



8210.002 

8220.002 

8230.002 

8240.002 

8250.002 



8260.002 

8560.001 

8570.001 

8580.001 

8590.001 



8630.001 

8690.001 

8700.001 

8730.001 

9350.001 



9350.002 

9410.001 

9530.001 

9590.001 


T 


710.001 

910.001 

1240.001 

2340.001 

2940.001 



3250.001 

3300.001 

5640.001 

5790.001 

5900.001 



5910.001 

6070.001 

6080.001 

6240.001 

6250.001 



6430.002 

6431.002 

7430.001 



TO 


6630.001 

7020.001 

7120.001 

7430.001 

8030.001 

Tl 


6640.001 

6900.001 

7120.001 

7430.001 


T2 


3690.001 

3710.001 

4140.001 

4220.001 

6810.001 



7020.001 

7420.001 

7870.001 

9080.001 

9110.001 



9140.001 

9170.001 

9180.001 



T3 


7060.001 

7160.001 




T4 


7110.001 

7120.001 
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T4 

7110. 0C1 

7120.001 




T5 

7120.001 

7160.001 

743U.001 

7940.001 


T6 

5620.001 

5640.001 

5790.001 

5900.001 

5910.001 


6070.001 

6080.001 

6240.001 

6250.001 

6430.002 


6431.002 

7060.001 

7420.001 

7940.001 


T7 

1430.002 

1430.003 

1590.001 

8280.001 


T8 

650.001 

660.001 

670.001 

670.002 

680.001 


690.001 

700.001 

710.001 

720.001 

730.001 


740.001 

750.001 

760.001 

770.001 

780.001 


790.001 

1080.001 

1340.001 

1580.001 

1800.001 


2300.001 

3390.001 

3410.001 

3430.001 

3440.001 


3450.001 

3460.001 

3470.001 

3480.001 

3490.001 


3500.001 

3510.001 

3520.001 

3530.001 

3540.001 


3550.001 

3560.001 

3570.001 

3580.001 

4060.001 


5650.001 

6030.001 

6200.001 

6431.001 

6502.001 


7330.001 

7420.001 

7430.001 

7450.001 

7490.001 


7600.001 

8150.001 

8170.001 

8180.001 

8190.001 


8200.001 

8210.001 

8220.001 

8230.001 

8240.001 


8250.001 

8260.001 




T9 

6380.001 

7430.001 

7940.001 

9700.001 


U 

I 80.001# 

4470.001 

4710.001 

4720.001 

4770.001 


4850.001 

4860.001 

4910.001 

5070.001 

5080.001 


5120.001 

5143.001 

5210.001 

5330.001 

5480.001 


5490.001 

9380.001 

6790.001 

7490.001 

7740.001 

9370.001 

UO 

9350.002 

9360.001 

9370.001 

9380.001 

9530.001 


9540.001 

9600.001 

9550.001 

9560.001 

9580.001 

9590.001 

U1 

9350.002 

9620.001 

9370.001 

9380.001 

9590.001 

9600.002 

U2 

3710.001 

3780.001 

7230.001 



U4 

4140.001 

4150.001 

4220.001 

4230.001 


U5 

4150.001 

8030.001 

4230.001 

6810.001 

7020.001 

7430.001 

U6 

4880.001 

4890.001 




U9 

4890.001 

4910.001 
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Detailed Program Description 

This section describes in more detail the program operations, keyed to the 
line numbers given above. Generalized descriptions are provided where a 
combination of the previous description and the comments embeded in the list- 
ing appear adequate; for the less obvious functions, more detail is provided 
here. To facilitate the description, it will be assumed that the program is 
memory resident, ready for execution. Program flow for Case #1 will be 
examined to illustrate operation. 

After the array declaration statements at line numbers 70 - 90, two output 
files are opened and margins axe defined for each (lines 130 - 160). File 
"OUTPUT. DAT” is used for 132 co 1 mm printer data, and "CRPLOT.MT" is an 80 
column file used for cross-plot data accumulation in the off-design point 
analysis. 

DATA statements 170 - 250 contain the 9X8 coefficient matrix elements 
used in the blown propeller aerodynamic coefficient lookup table subroutine. 

The matrix itself, J(9,8), is filled during the execution of statements 260 - 
300. The rows of J( ) represent angle of attack intervals of three degrees 
each. They start with the interval (-12,-9) degrees and end with the interval 
(12,14) degrees. The first four columns of J( ) represent linearized 
coefficients used in the lift coefficient equation and the second four columns 
represent the coefficients for the drag coefficient equation. The subroutines 
at 9200 - 9630 calculate the lift and drag coefficients as a function of angle 
of attack and momentum coefficient as 

C A » J(M) + J(*,2)C y + ^J(*,3) + J (V)C y ) (A14) 


C d «J(*,5) + J(*,6)C t ♦ ^J(*,7) + J(*,8)C y ^ a ’ a £^ (A1S) 

■ momentum coefficients, 

■ angle-of -attack interval lower limit 

■ angle-of-attack interval upper limit 

14 degrees, - 0.001 and C. ■ sin a. 
a a 
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The DATA statements between lines 310 - 500 are used to define the 
programmed cases 1-10. Each case requires thirty-one data elements. The 
data defined here is read into a 10 X 31 array Q(10,0:30) with statements 510 - 
550 in the order specified by the symbols list preceeding this program listing. 

Program lines 560 - 610 reset certain flags and accumulators specified in 
the comment statements. Lines 620 - 640 starts user interaction by asking for 
either a predefined case (1 - 10) or a signal for manual input (0). If a non- 
zero case is entered, airspeed, engine RPM, etc. are loaded from the Q( ) 
matrix at lines 680 - 780; if a zero is entered, the program interactively 
obtains the required data at lines 800 - 950. For the example Case #1, lines 
680 - 780 make V = 270, N1 -= 41.6, D = 6, X = 324, H2 = 10000, b = 3, PI = 0, 

C5 = 0, B7 = 1, C6 = 0, and C7 = 3. Between lines 950 - 1070, radian/degree 
conversion factors are defined and the atmospheric characteristics subroutine 
is invoked. Lines 1090 - 1200 interactivly query for additional case data; 
this information has already been supplied for the cases 1-10, so program flow 
jumps to line 1210 for the Case #1 example. 

The section from line 1210 to line 1330 calcuates run data based on the 
specific case under analysis. Variable descriptions are given in the REMark 
statements and equation references are also provided there. As mentioned 
earlier, equation references enclosed with angle brackets < > are from 
Reference 4, while those enclosed with square brackets [ 1 are from this 
report. Line 1340 branches around the case examination question at lines 1350 
- 1400, and lines 1420 - 1440 set up the control variables V4, T7, and V9 
according to which propellers are to be designed/evaluated. For Case ?/l, both 
non-blown and blown propellers will be designed (C7 =3); in this situation, V4 
* 1, T7 » 1, and V9 is left at its value of 0 defined at line 570. From the 
program variable list prior to the program itself; the control variables car. be 
interpreted as (1) the current case is non-blown (V4 = 1), (2) both non-blown 
and blown propellers will be evaluated (T7 - 1), and (3) output from the 
Analytic design section has not been completed (V9 = 0). When both non-blowr./ 
blown analysis is conducted in a single Case, the non-blown evaluation is 
always done first. 

Lines 1450 - 1570 inquire as to whether certain diagonistic output is to 
be included in the print file "OUTPUT. DAT'*, and whether during the induced 
velocity iterations (inflow calculations) certain data is to be routed to the 
user terminal (CRT). Responses to the two questions cause control variables V8 
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and W2 to be set to either a 1 (include output) or a 0 (no output). The 
selectable output can easily be located in the listing by first referring to 
the Symbol Cross-Reference Table just after the program listing. In this table 
the control variables V8 and W2 can be located, and all line numbers that 
include them can be found. 

Lines 1530 and 1590 route program flow around interactive entry for blade 
angle of attack distribution (lines 1600 - 1780) if either a defined entry lias 
been made, or if a manual entry non- blown/ blown case iu in progress and the 
first part (non-blown) has been completed (V4 = 0). For defined cases, lines 
1790 - 1810 load the A0A distribution from the Q( ) matrix. Lines 1820 - 2140 
set up the angle of attack intervals and the linearized coefficients for the 
lift and drag coefficient equations used in the non-blown propeller aerodynamic 
coefficient subroutine (located at lines 8410 - 8770). The equations realized 
in this subroutine are similar to those used In the blown propeller lookup 
table, and can readily be determined by inspection of the code. 

At line 2190, program flow is diverted to line 2540 for the blown 
propeller examination. For the example Case #1, the non-blown propeller is 
first designed so flow continues tc lines 2200 - 2240 where the non-blown 
propeller aerodynamic coefficients are obtained. Both analytic and strip 
integration design is effected for the non-blown propeller; if the analytic 
design has been completed (V9 = 1), line 2250 routes execution to line 2540 to 
avoid the first set of analytic printer output (lines 2270 - 2530). Between 
lines 2540 and 3120, the analytical design is made based on equations listed in 
the REMark statements. This design is carried out for both analytic end strip 
integration evaluation since some of the results obtained in the analytic 
calculations are used later by the strip integration section. Lines 3120 - 
3300 send more analytic output to the printer; control statements at lines 
3160 and 3260 route flow around PRINT statements if the analyti" output is 
complete . 

Lines 3370 - 3410 do the analytic chord and beta distribution calcuations 
for all cases 1, 2, 3, or 4. Lines 3430 - 3590 save or restore these chord and 
beta distributions for later case use. For the example Case #1, line 3430 
sa/es the beta distributions in the 5th column of the T( ) array if the 
non-blown propeller is being analyzed (V4 - 1). For Case #2 ( ins 3440), the 
chord distributions used are obtained from the 15th column of the T( array, 
which was saved at line 8180 during a Case #1 blown propeller design run. The 
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current scheme for saving and restoring data can be seen in the accompanying 
table; data is stored/read in the vicinity of lines 3430 - 3590 or 8160 - 8270. 
The rationale for this placement is that the 3400 locations occur before any 
strip integration calculations take place that are impacted by the chord or 
beta distributions, and that the 8100 locations occur after all strip integ- 
ration calculations have taken place and the data can be used by later cases. 


Case # 

Blown ? 

Action 03400’s 

Action 08100’s 

1 

no 

T(5,I)=Beta 

T(ll,I)=Chord 

1 

yes 

none 

T(15,I)=Chord, T(9,Il=Beta 

2 

yes 

Chords=T(15,I) 

T(12,I)=Chord, T(6,I)=Beta 

3 

no 

T(7,I)=Beta 

T(13,I)=Chord 

3 

yes 

none 

T(16,I)=Chord, T(10,I)=Beta 

4 

yes 

Chords=T(16,I) 

T(14,I)=Chord, T(8,I)=Beta 

5 

yes 

Beta=T(9, I) , Chords=T(15,I) 

none 

6 

yes 

Beta=T(6,I) , Chords=T(12,I) 

none 

7 

yes 

Beta=T( 10,1), Chords=T( 16 , I ) 

none 

8 

yes 

Beta=T(8,I) , Chords=T(14,I) 

r :>ne 

9 

no 

Beta=T(5,I) , Chords=T(ll,I) 

none 

10 

no 

Beta=T(7,I) , Chords=T(13.I) 

none 

Lines 

3600 - 3860 

consists of the rest of the 

analytic calculations 


vdescribcd in the REMark statements) and the last of the analytic print 
statements. Control statements at 3600 and 3800 route flow around these 
o.’wput statements if the analytic outpu for this case has already been done. 
Statement 3880 sets tne control variable V9 = 1 to indicate the analytic 
output has been completed for this example Case #1. 

Lines 3890 - 5620 implement the strip integration equations contained in 
the main part of this report or (for the induced velocity iterations) in this 
appendix. Details found in the REMark statements should prove ample to follow 
Che program flow through this section. Two areas will be expanded here for 

clarity. First, at line 4060 is the mechanism to have either a blown or a 

non-blown propeller. If the case identification is 1, 2, 3, or 4, and if the 
propeller under evaluation is blown (V4 = 0), then the jet velocity at the tip 

(station 19) is set to 95% of the speed of sound; otherwise, the jet velocity 

is set to and “lie propeller is non-blown. 

Secondly, if the propeller is blown, it is blown from root to tip (station 
1 to station 19). If partial span blowing is desired, then coding changes in 
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this section are necessary to (1) reduce the last blade station number at which 
blowing will occur, and (2) match the lift coefficient at the transition 
between the blown and non-blown blade sections. The following code replacement 
could be used as a basis for a partial span evaluation, with the Liner 2/3 
(root to station 14) blown, and the outer 1/3 (station 15 to 20) non-blown. 

Replacement Code for Evaluation of Partial Span Blowing 


3960 

3970 

3980 

3990 

4000 

4010 

402C 


****** 


REM** ** 

REM** ROUTE PROGRAM FLOW BASED ON BLOWN OR NON-BLOWN OPTION ** 

REM** (BLOWN, V4=0; NON-BLOWN, V4=l) ** 

IF V4=l THEN GO TO 4670 

REM** ** 


4030 

4040 

4050 

4060 

4070 

4080 

4090 

4100 

4110 

4120 

4130 

4140 

4150 

4160 

4170 

4180 

4190 

4200 

4210 

4220 

4230 

4240 

4250 

4260 

4270 

4280 

4290 

4300 

4310 

4320 

4330 

4340 

4350 

4360 


REM********** DETERMINE MOMENTUM COEFFICIENTS AT STATION 14 ************ 
REM********** NECESSARY TO MATCH LIFT OF S/C AIRFOIL ************ 
REM 
1=14 

G0SUB 9470 

U3=(C(1,14)-U(U0,1)-U1*U(U0,3))/(U(1?0,2)+U1*U(U0,4)) 

REM 

REM********** U3 IS THE MOMENTUM COEFFICIENT AT STATION 14 ************ 
REM 

REM********** NOW CALCULATE THE JET VELOCITY AT STATION 14 **** 

REM 

H( 1 , 14 ) =V( 14 ) *26 . 4*U3** . 5 
REM 

REM********** THE REQUIRED JET PRESSURE AT STATION 14 IS 
REM 

H( 2 , 14)=0 . 5*H*H( 1 , 14)**2 . +P2 
REM 

REM********** THE HUB PRESSURE TO GIVE THIS STATION 14 IS 
REM 

U4=( ( PI*N1*K( 14 ) *D ) **2 ) / ( 3432*T2 ) 

U5=H(2,14)/EXP(U4) 

REM 

REM******** ’ IF REQUIRED HUB PRESSURE IS LESS THAN STATIC ************ 
REM*******- PRESSURE, USE STATIC PRESSURE AT HUB! ************ 

REM 

IF (U5CP2) AND (C6<fc) 

THEN PRINT "REQUIRED HUB PRESSURE IS LESS THAN STATIC FOR CASE #"}T8 
IF (U5<=P2) AND (C6<>1) THEN N1=N1*(T/T6)**. 5 
IF U5<P2 THEN U5=P2 
REM 

REM*********************************************** - *********** ********** 



REM* CALCULATE JET PRESSURE AT EACH BLADE 

REM 

OR 1=1 TO 14 


STATION [EQN. 6] 

***• V*>f .';************* 
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4370 U4<(PI*N1*K(I)*D)**2)/(3432*T2) OF POOR QUALITY 

4380 H( 2 , 1 )=U5*EXP(U4 ) 

4390 REM ^ — a; 

4410 REM* CALCULATE JET VELOCITY AT EACH STATION ALSO [EQN. 2] 

4430 REM 

4440 H(1,I)=(((H(2,I)-P2)*2)/H)**.5 
4450 REM 

4460 R ]^**** * ****************************************** v "”^^^ 

4470 REM* * 

4480 REM* CALCULATE EACH BLADE STATIONS LOCAL VELOCITY [EQN. 7] * 

4490 RcM* * 

4496 

4497 H( 3 , I ) = ( V**2+( PI*K( I ) *D*N1 )**2 ) ** . 5 

4498 V( I )=SQR(V**2+( 2*PI*K( I )*R6*N1 )**2) 

4499 REM 

4500 REM ***** * ****************************** ***VnV*VnVVr;V***Vf**Vrfr.'r**-A-**-. .■******* 

4501 REM* ... AND FINALLY, CALCULATE THE MOMENTUM COEFFICIENT 

4502 REM* AT EACH BLADE STATION [EQN. 1] 

4503 REM ******************************************************************** 

4504 REM 

4505 REM 

4506 U(l,I)=((H(l,I)/H(3,I))**2)/696.2 

4507 REM**********' ntnVVr ' v *******' v *** v ^ 

4508 REM* END OF LOOP 

4509 REM 

4510 NEXT I 

Other changes may be required to the program depending on the method used 
to match developed thrust to required thrust. For example, if the momentum 
coefficients are scaled -;o achieve the required thrust, then program flow after 
jet adjustment must return to statement 4120 to recalculate the required hub 
pressure and other affected quantities. In addition, the blown propeller 
aerodynamic coefficient subroutine must be modified at line 9330 to account for 
the partial span blowing to station 14. Calls to the non-blown propeller aero 
coefficient routine must set V3 to station 15, the first non-blown location. 

Lines 5630 - 6510 test the developed thrust against the required thrust, 
and if they do not match within 1% of the required thrust, some form of 
adjustment is implemented. The specific form of adjustment employed was 
discussed earlier, in the Case Table Identification of Table Al. For this 
program, Alpha, Beta, engine speed, or chord scaling is used to obtain the 
required values of thrust, with the specific parameter controlled by the values 
of 87, C5, C6, or Case #. For the non-blown propeller of example Case #1, the 
chords are scaled at line number 6430, because (1) B7 = 1 causes program flow 







to skip the Beta adjustment (line 5860), (2) V4 = 1 causes a skip of the Alpha 
adjustment (line 6030), (3) Case # 1 (T8 = 1) causes a skip of the engine speed 
adjustment (line 6200), and (4) C5 <> 1 and Case # < 5 results in the execution 
of line 6430, which is the chord adjustment. Similar logic can be applied to 
verify tha form employed by other cases. After any adjustment is made, program 
flow returns to line 4560 to recalculate the induced velocity component. This 
iterative process occurs until either a satisfactory thrust match has been 
obtained, or the user terminates the case evaluation. 

Lines 6520 - 7290 calculate torque, air mass flow, compressor power 
requirements, etc. This section has ample comments to follow the operations 
with references as appropriate. Line 7270 calls a subroutine to determine 

the horsepower available at a^specific engine speed and full throttle for a 
Turbocharged engine of 520 in' 3 displacement. 

Lines 7300 - 7520 output cross-plot data to the " CRPL0T.DAT" file for the 
off-design cases (5 - 10). Each time a thrust adjustment is made in the 
off-design runs, the cross- plot data file has data added to it. In this 
fashion, data is accumulated as a function of the parameter making the thrust 
adjustment. This data is then used to generate plots like Figure 4 in the main 
report. The variables seen in lines 7420, 7430, 7450, and 7490 are saved at 
each data point; the interpretation of these variables can be made by referring 
to the List of Symbols located just before the program listing. There are four 
record types output to this file; the record type identification is always the 
first two-digit number in each record. Valid record types are 00, 01, 02, 03, 
and 99. Record type 99 is a data point delimiter. Data analysis of the cross- 
plot file requires a knowledge of the output format; the easiest way to obtain 
that information is to run the program and then type (part of) the cross-plot 
file. When data is added to this file, no printer file output is generated. 
This action is controlled by the variable A4, as can be seen in lines 7330 and 
7520. 

Lines 7530 - 8100 output printer data for the strip integration equations. 
The printer file is currently set up for a 132 column printer and essentially 
all columns are required for the blown propeller output quantities. Fewer 
are needed for the non-blown propeller output, but an 80 column printer is 
still not adequate. Control variable VI is used to select the correct output 
statements (non-blown/blown). 

Lines 8110 - 8270 save final chord and beta values as was discussed 
earlier. Lines 8280 - 8300 control the program flow for the multiple 
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evaluation cases, such as the example Case # 1. This case first requires the 
analytic design of a non-blown propeller, followed by the strip integration 
design of the same non-blown propeller. After a % output is complete (line 
8100), and the chord values have been saved, statement 8290 sets control 
variable V4 (which was 1 for the non-blown propeller) to a 0 (for a blown 
propeller design). Control then returns to statement 1790 to design/evaluate 
the blown propeller. When execution returns to statement 8280, V4 = 0 causes 
the program flow to resume at line 8320. 

Lines 8320 - 8400 execute at the end of each case evaluation and allow 
the user to run another case or to stop the program. If another case is 
selected, control returns to line 570; otherwise the program stops at line 
8400. 

Most of the subroutines found between lines 8410 and 9710 have already 
been discussed. The one exception to this is the Simpson’s Rule numerical 
integration subroutine found at 8780 - 8920. This routine has the integrand 
passed to it in the Z( ) vector which has 20 elements. Lines 7480 to 7500 
sum the odd numbered Z( ) elements into the temporary variable Q, that is 
Q = Z(l)+Z(3)+Z(5)+Z(7)+Z(9)+Z(ll)+Z(13)+Z(15)+Z(17)+Z(19) 
and lines 7520 - 7540 sum the even numbered Z( ) elements into the temporary 
variable R, as 

R = Z(2)+Z(4)+Z(6)+Z(8)+Z(10)+Z(12)+Z(14)+Z(16)+Z(18)+Z(20) 

These quantities are then formed into the integral by application of 
Z = 0.05 (4Q + 2R) / 3 (where Z is a scalar) 

The integral is returned to the calling program in the scalar variable Z. 

Suggestions for Tailoring the Program 

To start a new effort using this program requires modification of two main 
program elements: (1) the airfoil and engine characteristics, and (2) the pre- 
defined case identification, and program logic flow. The case identification 
and program flow have been discussed earlier. The remaining changes are 
discussed here. 

The main task in adapting this program to use some arbitrary airfoil shape 
and engine characteristic is to obtain the piecewise-linear coefficients used 
by the lookup table subroutines (Lines 8410 - 8770, 9200 - 9630, and 9640 - 
9710). This data is obtained from airfoil lift/drag polars by first segmenting 
the lift/drag curve into three degree parts (for the current program), and then 
obtaining coefficients for the slope-intercept form of a straight line best 
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representing the lift/drag curve over each angle of attack range. This data is 
then stored in the DATA statements (170 - 250) used to load the J( ) array for 
the blown propeller, or in statements 1990 - 2140 for the non-blown propeller. 
Engine characteristics are derived in a similar manner, and the slope-intercept 
data used in statement 9700 for the horsepower available subroutine. 

The remaining comments are intended to assist in the conversion of this 
program to run on some other machine. Since it is impossible to consider all 
implementations of BASIC that might be encountered, the statements that appear 
to be non-standard are discussed. Establishing new cases for evaluation, 
modification of the thrust matching procedure to account for other parametric 
variation, reducing the blowing span, and changing the engine characteristics 
have all been discussed. 

Lines 70 - 90 will cause problems with those versions of BASIC that do not 
allow multiply subscripted arrays. This poses a very difficult problem, and 
most likely, makes the program unusable, fortunately, not many versions have 
such limitations. Lines 130 - 160 may cause problems in the form used to open 
files. A printer disk file as such may not be required} if not, then replace 
line 130 with the appropriate type. The MARGIN statement may not be allowed; 
delete it, or find a suitable replacement to tell«the computer how wide the 
longest line will be in each file. 

Line 620 prints the current date on the user's terminal; it can easily be 
removed if the function is not otherwise available. The majority of the 
program is "standard" BASIC, so the next possible problem area are in the 
output sections scattered throughout the program. To a large extent, the PRINT 
USING statement was utilized to more precisely control the output format. If 
the selected version of BASIC has a PRINT USING statement, it may not agree 
exactly with the one used in the current program, but only small changes should 
be required. If no PRINT USING is available, then all strip integration output 
will have to be redone to fit whatever PRINT statements are available. The 
format variable names used (e.g. F0R1$) may also be a problem; if so, change 
the names. String concatenation is used co achieve long format strings (longer 
than 30 characters) as can be seen in lines 7360 and 7400. If this cannot be 
accomplished in this manner, some other scheme must be devised. 

It is likely that the program size will cause difficulties in transferring 
this program to another machine. Since the machine used for the current study 
compiled the BASIC source program to its own machine code, and used common 
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memory resident run-time libraries, the actual number of bytes of memory 
required for interpretative machine execution is not known. Obvious techniques 
to reduce memory requirements include (1) removal of the REMark statements and 
the comments on each line, and (2) segmentation or chain operation of the 
program. In an interpretative BASIC computer, REMark statements require 
available user memory just as executable statements do, but they do not affect 
program results. If segmentation or chain operation is used, the program 
can be divided into parts in a fashion similar to the way the program was 
discussed in the General Description portion of this appendix. 

Concluding Remarks 

A comprehensive discussion of the propeller design program used in the 
current study has been given. This discussion first addressed the general 
characteristics of each major program section, and then described in detail the 
manner in which predefined cases were set up. A complete program variable 
list, with the engineering units used by the program (as appropriate), was 
given just prior to a complete program listing. A symbol table (program 
variable) cross-reference listing was then given to facilitate program 
understanding. A detailed program walk-through followed these listings, keyed 
to the program line numbers, and using an example case for clarity. The 
appendix ends with some comments on adapting the program to some arbritary 
airfoil and engine characteristic and transferring the program to other 
machines . 
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